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ABSTRACT 

The  purpose  of  this  thesis  is  to  investigate  the  state 
of  the  art  of  radar  signal  design  as  well  as  radar  signal 
processors  and  determine  the  actual  trends  in  modern  radar 
design.   The  use  of  a  digital  general  purpose  radar  signal 
processor  is  discussed.   The  concepts  of  ambiguity  and  auto- 
correlation function  are  investigated  in  regard  to  radar 
resolution  capabilities.   The  concept  and  analytical  devel- 
opment of  the  DFT/FFT  are  presented.   Quantization  noise  in 
a  digital  MTI  processor  and  its  effects  in  the  improvement 
factor  are  analyzed.   Optimization  techniques  for  the  response 
curve  of  digital  MTI  processors  using  staggered  PRF  are  in- 
vestigated.  The  SAR  concept  and  analysis  as  well  as  techniques 
to  obtain  low  correlator  rates  in  the  SAR  digital  processors 
are  presented. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 8 

II.  RADAR.   HISTORY  AND  APPLICATIONS 10 

III.  SIGNAL  PROCESSING  IN  RADAR 12 

A.  INTRODUCTION  TO  THE  RADAR  RANGE  EQUATION 12 

B.  RELATION  BETWEEN  THE  RADAR  RANGE  EQUATION 

AND  SIGNAL  PROCESSING 14 

C.  USE  OF  COMPUTERS  AS  SIGNAL  PROCESSORS  IN  RADAR---  15 

IV.  VARIOUS  ASPECTS  OF  THE  THEORY  OF  RADAR 

SIGNAL  PROCESSING 20 

A.  SIGNAL  PROCESSING 20 

B.  AMBIGUITY  FUNCTION  AND  AUTOCORRELATION 23 

1.  Concept 23 

2.  Range  Ambiguity  Function 23 

3.  The  Velocity  Ambiguity  Function 26 

C.  PULSE  COMPRESSION 30 

1.  Concept 30 

2.  Linear  FM  Chirp  Pulse  Compression 51 

3.  Matched- Filter  Approach 36 

4.  Use  of  Discrete  Frequency  Sequences 

in  Pulse  Compression 37 

D.  MOVING  TARGET  INDICATOR  (MTI) 45 

E.  DIGITAL  SIGNAL  PROCESSING,  DFT  AND  FFT 59 

F.  SYNTHETIC  APERTURE  RADAR 73 

V.  RECENT  DEVELOPMENTS  IN  TWO  MAJOR  AREAS 86 

A.  DIGITAL  MTI 86 

B.  DIGITAL  SAR -109 


VI.   CONCLUSIONS.  TRENDS 123 

BIBLIOGRAPHY 124 

INITIAL  DISTRIBUTION  LIST 126 


ACKNOWLEDGEMENT 

The  author  wishes  to  acknowledge  the  support  and 
encouragement  provided  him  by  Associate  Professor  John  M. 
Bouldry.   In  addition,  the  author  wishes  to  thank  all  the 
Dudley  Knox  Library  staff. 


I.   INTRODUCTION 

The  processing  of  signals  has  been  always  a  necessity 
whenever  information  was  to  be  transmitted  through  a  channel 
between  a  source  and  a  destination.   It  was  not  until  1948 
that  Shannon  presented  the  basic  theory  of  information  re- 
lating the  source  entropy  with  channel  capacity  and  the 
probability  of  error.  Techniques  of  signal  processing  were 
already  in  use,  but  the  results  given  by  the  information 
theory  clarified  the  limits  that  could  be  achieved. 

Early  analog  processors  with  enough  bandwidth  and  reason- 
able signal  to  noise  ratio  could  process  data  at  extremely 
high  data  rates,  but  the  components'  low  stability  and  high 
cost,  and  sometimes  low  versatility,  were  responsible  for  a 
gradual  substitution  of  analog  by  digital  processors  when 
the  price  of  digital  logic  went  down,  and  the  digital  hard- 
ware technology  suffered  a  big  jump  in  the  last  ten  years. 

In  fact,  in  today's  radar  signal  processing  the  trend  to 
a  substitution  of  analog  by  digital  processors  is  inevitable. 
The  theory  and  boundaries  of  digital  signal  processing  are 
continually  being  improved  as  well  as  hardware  components. 
Speed  is  almost  no  longer  a  problem  and  versatility  is  immense. 

Chapter  II  is  an  introduction  to  radar  history.   In  Chapter 
III  signal  processing  is  related  with  the  radar  problem  through 
the  radar  range  equation  and  a  general  digital  processor  for 
a  radar  system  is  presented.   In  Chapter  IV  an  analysis  is 
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made  of  the  present  signal  processing  techniques  in  radar 
with  special  emphasis  into  digital  approaches.   In  Chapter  V 
an  investigation  is  made  into  the  recent  developments  and 
problems  of  two  major  areas  of  radar  digital  signal  process- 
ing.  In  Chapter  V  an  overall  view  of  the  achievements  and 
trends  of  radar  signal  processing  is  explored. 


II.   RADAR.  HISTORY  AND  APPLICATIONS 

The  electromagnet  ism  and  electromagnetic  wave  propagation 
theories  of  Maxwell,  later  (1886)  experimentally  verified  by 
Hertz,  contain  all  the  necessary  background  to  understand  the 
principles  of  radar.   The  first  experiments  on  detection  with 
radio  waves  started  in  1903.   Due  to  inadequate  technology 
the  results  were  very  poor.   The  obtained  ranges  were  less 
than  those  achieved  with  optical  systems  and  this  was  enough 
to  discourage  any  interest  in  pursuing  the  experiments.   It 
Avas  only  in  1922  at  the  Naval  Research  Laboratory  that  a 
wooden  ship  was  detected  using  a  CW  radar  with  the  transmitter 
and  the  receiver  as  separate  units.   From  then  on  an  increas- 
ing interest  in  radar  technology  became  evident.   In  1930 
using  a  bistatic  radar,  the  first  aircraft  was  detected,  and 
by  1932  the  detection  ranges  were  already  in  the  50-mile 
region.   By  that  time  the  pulse  techniques  were  not  yet 
explored,  so  the  information  was  in  the  presence  or  absence 
of  a  target  and  not  in  range  information.   But  in  1935  the 
British  successfully  used  pulse  techniques  to  measure  distance 
to  the  targets.   In  1936  the  detection  range  was  already  in 
the  90-miles  region.   The  basic  principles  were  understood. 
Its  next  steps  were  in  technology  improvements  in  order  to 
get  stable  elements  and  more  power  output  to  increase  the 
range  coverage.   Higher  frequencies  were  used  in  order  to 
work  with  smaller  elements  and  high  gain  antennas. 
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The  discovery  of  the  magnetron  was  an  important  step  since 
powers  went  up  by  a  factor  of  100,  and  wavelengths  of  10  cm 
were  obtained.   But  the  greatest  incentive  for  further  devel- 
opments in  the  radar  field  came  from  the  military  necessities 
of  World  War  II.   After  the  War  a  stagnation  of  about  five 
years  slowed  down  the  rhythm  of  improvements  in  the  field. 
But  in  the  fifties  the  introduction  of  the  high  power  klystron 
increased  not  only  the  power  output  but  the  frequency  stability 
necessary  for  coherent  detection  as  well  as  coherent  MTI ;  also 
very  low  noise  receivers  were  implemented. 

From  then  on,  with  more  and  more  perfect  technology,  more 
complex  systems  appeared.   The  introduction  of  the  computer 
as  a  storage  and  control  element  was  of  extreme  importance. 
All  ballistic  systems  and  very  large  coverage  systems  suffered 
a  big  impulse.   Digital  systems  improved  the  implementation 
of  stereable  array  antennas,  synthetic  aperture  radar,  track- 
ing systems,  as  well  as  all  digital  processing  techniques. 

With  the  development  of  radar,  the  statistical  nature  of 
the  detection  problem  gave  rise  to  an  increasing  interest  in 
the  study  of  the  statistical  properties  of  clutter,  radar 
cross  section,  rain,  etc.   So,  due  to  a  more  precise  knowledge 
(modeling)  of  the  systems  environment,  the  radar  systems 
changed  depending  on  the  type  of  application. 

Thus,  a  design  to  optimize  a  given  application  seldom 
can  be  used  in  others.   Today  the  radar  design  field  is  wider. 
Applications  vary  from  military  to  civilian  to  scientific. 
With  improvements  in  technology  and  more  accurate  modeling, 
research  in  the  field  is  far  from  becoming  saturated. 
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III.   SIGNAL  PROCESSING  IN  RADAR 

A.   INTRODUCTION  TO  THE  RADAR  RANGE  EQUATION 
Only  the  analysis  of  the  equation 

R  =  £(ni,n2,  ...  i^)  (3-1) 

where  R  is  the  distance  from  the  radar  antenna  to  the  target, 
and  (n, ,n~,  ...  nn)  are  the  radar  and  environmental  parameters, 
can  give  the  necessary  knowledge  to  maximize  R.   If  P   is  the 
peak  pulse  power  radiated  and  G  the  gain  of  the  antenna,  the 
transmitted  power  density  is 

.   PtG 


4  7T  R2 


and  the  received  power  density,  after  reflection  in  a  target 
of  cross  section  a   will  be 

P  G  a 

(3-2) 


(4ttR2)2 

If  (3-2)  is  multiplied  by  the  antenna  aperture  A  ,  the 

received  power  will  be 

P.  G  a.  A 
p   =  -1 L_£  (3_3) 

r    (4ttR2)2 

If  S  .   is  the  minimum  received  signal,  then  the  maximum 

distance  is 

P.  G  AQ  a  + 
it       t    e  t 


(4tt)  S 


r^    =  _i E — ^_  (3-4) 

5min 
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Since  by  definition  of  noise  figure  (F  )  of  a  receiver 


S.  =  N.  F   (  V^) 

1     1   n  v  N0  ' 


i  =  input 
o  =  output 


where 


Ni  ■  kTBn 


input  noise 


Then,  combining  (3-4)  with  (3-5)  and  (3-6) 


R4 
max 


PtGAe°t 


(4Tr)2kTB    F    fSn/Nn)     - 
*•      >  n   nv    o'    o^min 


(3-5) 


(3-6) 


Since  in  a  pulse  radar  the  average  power  is 


P   =  Pfxf 
av     t    r 


x  =  pulse  width 
fr  =  pulse  rejection  frequency 


and  the  receiver  bandwidth  B   is  approximately  1/t 


R4 

max 


P       G  A    a  + 
av        e     t 


(4i02kTf    F    (S0/N    )     . 
rnuo-/min 


(3-7) 


where  (S  /N  )     is  the  minimum  signal  to  noise  ratio  at  the 
v  o   o  mm  & 

output  of  the  linear  section  of  the  receiver,  necessary  for 
detection.   The  first  important  characteristic  in  equation 
(3-7)  is  the  fact  that  some  of  the  variables  can  be  considered 
deterministic  but  others  must  be  analyzed  as  random  variables. 
The  random  variables  are  a.    and  S/N.   Statistical  discriptors 
for  a   are  already  extensively  studied  [1-2]  not  only  in  the 
general  case  but  also  for  some  specific  cases,  as  for  different 
types  of  clutter.   Since  the  noise  is  normally  treated  as 
Guassian,  all  necessary  statistical  discriptors  are  available. 
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The  statistical  discriptors  together  with  the  decision 
criteria  will  determine  the  probability  of  detection  and 
probability  of  false  alarm. 

B.   RELATION  BETWEEN  THE  RADAR  RANGE  EQUATION 
AND  SIGNAL  PROCESSING 

After  the  War  (1945)  most  of  the  effort  in  increasing 

range  capabilities  (eq.  #3-7)  was  focused  on  increasing 

average  power  (P   ) ,  antenna  gains  (G) ,  and  mixers  and 

receivers  with  lower  noise  figures  (F  ) .   With  the  magnetron 

and  the  klystron,  the  peak  powers  reached  the  megawatts  range. 

Soon  it  was  verified  that  since  the  technology  was  already  in 

the  limits  of  transmitted  peak  power,  any  increase  in  peak 

power  was  questionable  from  a  technical  and  a  financial  aspect, 

Also,  since 

1/4 
R  =  kP  ,  k  =  constant 

In  R  =  In  k  +  j   In  P 

d_R  =  1  dP 
R    4  P 

That  means,  that  the  percent  gain  in  range  is  only  1/4  of  the 
percent  gain  in  power.   The  pulse  width  could  not  be  increased 
much  more  because  of  the  required  range  resolution.   This  led 
to  the  development  of  what  today  is  called  signal  processing 
techniques.   Waveforms  were  studied  in  order  to  increase  the 
average  power  but  maintaining  range  resolution.   The  basic 
techniques  will  be  presented  in  Chapter  IV. 

The  brute  force  approach  to  increase  the  range  of  detec- 
tion (3-7)   was  used  until  1955.   The  technology  then  shifted 
to  an  improvement  of  the  signal  to  noise  ratio.   If  the  new 
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signal  to  noise  ratio  is 

(S/N)  =  D(S0/NQ),  D  >  1. 

Since  D  >  1,  there  is  an  increase  in  the  maximu...  range  of 
detection  (3-7) . 

C.   USE  OF  COMPUTERS  AS  SIGNAL  PROCESSORS  IN  RADAR 

As  range  and  aximuth  resolution  became  finer,  as  multi- 
target  processing  due  to  digital  techniques  became  more 
effective,  the  conventional  final  processor  in  the  radar 
system,  the  man,  gradually  was  substituted  by  the  computer. 
The  total  amount  of  information  and  the  necessary  speed  for 
processing,  tax  the  capabilities  of  man.   So,  digital  tech- 
niques are  gradually  being  substituted  for  not  only  those 
processors  that  used  analog  devices,  but  also  where  human 
operators  were  used. 

Figure  1  gives  an  overview  of  the  techniques  presently 
used  in  digital  signal  processing.   As  in  the  analog  theory 
linear  approximation  is  used.   So,  using  the  theory  of  linear 
discrete  time  invariant  systems  there  are  two  common  approaches 
digital  filtering  and  spectral  analysis.   Using  the  theory  of 
digital  filtering, algorithms  can  be  implemented  that  correspond 
to  finite  impulse  response  filters  (FIR)  or  infinite  impulse 
response  filters  (IIF).   Using  spectral  analysis,  two  methods 
are  possible:  implementation  of  Fast  Fourier  Transform  (FFT) 
algorithms  or  use  of  statistical  spectrum  analysis. 

Both  digital  filters  and  spectral  analysis  are  affected 
by  quantization  noise,  which  must  be  taken  into  account. 
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As  seen  in  Fig.  1,  the  radar  appears  as  one  of  the  prime 
applications  for  digital  systems. 

It  was  in  the  1960 ' s  that  the  first  low  data  rate  digital 
processors  were  implemented  in  radar  systems,  basically  to 
perform  tracking  and  weapon  order  computations.   The  next 
step  was  the  implementation  of  the  MTI  processor.   A  flexible 
programmable  digital  processor  was  not  achieved  until  very 
recently.   The  low  cost  of  semi-conductor  memory  and  the 
advanced  technology  in  microprogramming,  with  inherent  increase 
in  speed,  were  the  main  reasons  for  this  achievement.   Now  a 
typical  general  purpose  signal  processor  for  radar  will  be 
described.   Figure  2  shows  a  block  diagram  of  a  digital  signal 
processing  system  for  a  radar.   The  main  processing  unit  is 
the  arithmetic  pipeline,  which  is  a  general  purpose  signal 
processor  (GPSP) .   This  consists  of  five  sections  selected 
in  a  way  that  the  most  frequently  necessary  operations  can 
be  performed;  Fast  Fourier  Transform  (FFT)  ,  recursive  and 
nonrecursive  digital  filtering  such  as  MTI,  threshold  gener- 
ation, peak  detection  noncoherent  video  integration  and  range 
and  angle  estimation.   The  first  block  in  the  pipeline,  the 
matrix  switch,  selects  data  from  data  memories,  also  from  the 
receiver  through  a  track  buffer,  or  from  anoth  er  GPSP  via  a 
wrap-around  connection.   The  second  block,  data  scaling, 
connects  the  data  for  floating  point  arithmetic.   The  complex 
multiplication  block  does  the  equivalent  for  four  real  multi- 
plications and  two  real  additions,  which  corresponds  to  a 
complex  multiplication.   Since  some  operations,  like  integration, 
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The  re-normalization  block  is  sometimes  needed.   Connected 
to  the  arithmetic  pipeline  is  a  memory  section  and  a  control 
section.   The  output  of  the  arithmetic  pipeline  is  relayed 
to  an  input/output  control  unit  (IOC)  which  is  a  buffer  for 
the  radar  data  as  well  as  a  macro  control  for  the  processing 
functions  of  the  GPSP.   The  advantages  of  this  processing  system 
are  that:  once  the  data  is  started  the  process  is  continuous, 
and  does  not  follow  the  normal  sequence  of  fetch,  operate,  and 
store,  characteristic  of  general  purpose  computers;  the  pro- 
cessor has  a  horizontal  structure  as  opposed  to  the  normal 
vertical  structure  of  computers;  this  means  that  an  instruc- 
tion, after  execution,  causes  all  operations  specified  to  be 
executed  in  a  clock  cycle;  there  is  a  hardware  separation  of 
instructions  and  data  in  the  memory,  which  simplifies  the  job 
of  the  programmer;  the  control  memory  is  made  of  RAM's  which 
implies  a  much  greater  flexibility  in  subroutine  changes;  and 
the  execution  of  instructions  is  within  the  same  time  frame 
that  corresponds  to  a  range  or  doppler  cell. 

Four  memories  interact  with  the  arithmetic  pipeline.   Two 
of  them  store  and  recirculate  data,  the  third  stores  weighting 
coefficients  for  such  filters  as  MTI ,  FFT  or  pulse  compression, 
the  fourth  stores  micro  opcodes.   Data  sequences  and  arithmetic 
operations  are  both  controlled  through  software,  by  two  control 
units.   The  program  that  resides  in  the  sequence  control  memory 
(SEQ-CONTMEM) ,  informs  the  machine  of  range  and  doppler  dimen- 
sions of  the  problem.   The  program  in  the  Arithmetic  Control 
Memory,  ACM,  informs  the  machine  of  the  particular  algorithm 
to  be  used. 
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IV.   VARIOUS  ASPECTS  OF  THE  THEORY 
OF  RADAR  SIGNAL  PROCESSING 

A.   SIGNAL  PROCESSING 

The  desire  to  transmit  information  is  closely  related  to 
the  processing  of  signals.   Since  normally  the  information  to 
be  transmitted  cannot  flow  through  the  available  channels,  it 
is  necessary  that  it  be  processed  to  use  the  chosen  channel. 
Looking  at  Fig.  3,  the  inverse  process  has  to  be  accomplished 
in  order  to  get  the  information  in  an  understandable  form  at 
the  receiving  end.   The  process  is  not  so  simple,  since  in 
general  the  processors  themselves  and  the  channel  introduce 
noise  that  may  or  may  not  eliminate  or  change  the  understanding 
of  the  message.   So  there  is  a  need  to  interpret  the  results 
with  some  decision  criteria.   Here  decision  theory  plays  an 
important  role.   Even  for  the  same  type  of  application,  for 
instance,  radar,  decision  criteria  are  not  the  same.   The 
criteria  may  be  of  constant  false  alarm  rate,  or  may  be  of 
fixed  probability  of  detection.   Basically  every  application 
or  system  must  have  its  decision  criteria. 

There  are  two  main  approaches  to  the  processing  of  signals, 
the  time  domain  approach  and  the  frequency  domain  approach. 
Besides  the  fact  that  they  use  completely  different  hardware, 
in  the  pure  mathematical  way,  they  are  duals  of  each  other 
through  the  Fourier  transform  pair.   So  the  only  reason  for 
one  type  of  approach  is  due  to  more  perfect  reliable  type  of 
hardware . 
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With  the  advance  of  computer  techniques ,  both  approaches 
are  split  between  digital  and  analog  signal  processing. 

It  was  only  after  World  War  II  that  electronic  engineers 
became  interested  in  the  applicability  of  digital  hardware 
techniques  in  signal  processing  areas.   But  it  was  not  before 
the  late  sixties  that  digital  techniques  suffered  a  big  jump 
replacing  some  of  the  analog  techniques.   The  foundations  of 
digital  signal  processing  can  be  related  to  the  Laplace's 
Z- transform  theory;  but  only  in  the  mid  sixties  was  the  theory 
presented  in  a  formal  way.   Various  papers  appeared  at  that 
time  related  to  this  subject,  but  only  in  1969  was  the  first 
attempt  made  for  a  comprehensive  theory  of  digital  signal 
processing  [21],   Most  recently  (1975),  two  very  comprehensive 
books  by  Oppenheim  [3]  and  Rabiner  [4)  can  be  considered  as 
giving  an  excellent  treatment  of  the  subject,  not  only  through 
its  mathematical  structure  but  with  very  good  applications. 
Today's  trend  is  definitely  towards  the  substitution  of  analog 
processing  by  digital  processing.   The  speed  achieved  in 
today's  digital  processors,  the  very  fast  algorithms  used  to 
implement  the  FFT,  the  gradual  substitution  of  infinite  impulse 
response  methods  to  finite  impulse  methods  due  to  a  better 
knowledge  of  theory  and  higher  efficiency  of  calculations  are 
the  major  factors  responsible  for  the  increasing  shift  to 
digital  techniques.   There  are  yet  some  type  of  applications 
where  the  speeds  are  so  high,  or  the  digital  hardware  is  so 
complicated  that  analog  techniques  are  still  used. 
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B.   AMBIGUITY  FUNCTION  AND  AUTOCORRELATION 

1.  Concept 

Considering  only  non-accelerating  targets,  the  two 
kinds  of  information  of  interest  are  the  position  of  the 
target  and  its  velocity.   Given  two  targets,  the  capacity 
to  differentiate  between  their  positions  and  velocities  is 
vital  for  most  radar  applications. 

The  capacity  to  differentiate  in  range  is  called 
range  resolution,  and  in  velocity,  velocity  resolution. 

As  far  as  range  resolution  is  concerned,  it  is 
obvious  the  shorter  the  pulse  the  better  is  the  differentia- 
tion between  the  two  targets,  that  is,  the  higher  is  the 
range  resolution.   With  the  measurement  of  velocity,  since 
it  is  directly  related  to  the  phase  difference  between  sig- 
nals, the  longer  the  pulse  the  greater  the  number  of  cycles 
that  can  be  compared  and  consequently  the  better  the  doppler 
and  the  velocity  resolution.   So,  at  first  glance  it  looks 
like  to  optimize  one  of  the  resolutions  the  other  will  have 
to  be  sacrificed.   The  other  solution  is  to  try  to  get  an 
optimum  solution  for  both  cases  at  the  same  time.   There  is 
a  quantitative  way  to  express  these  conditions  in  a  precise 
way.   That  was  the  reason  for  the  appearance  of  the  mathe- 
matical concept  of  ambiguity  function. 

2 .  Range  Ambiguity  Function 

Let  the  transmitted  signal  be  represented  by  S(t)  = 
Re[\jj(t)],  i|'(t)  =  u(t)  exp.  j  oo  t.   Then,  with  no  doppler 
present  the  received  signals  from  two  different  targets  are 
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Y(t)  and  lp(t-x) 
t=  equivalent  time  difference  between  the  two  targets. 

If  it  is  chosen,  as  a  measure  of  the  resolution,  the  mean  square 
of  the  difference  between  the  two  received  signals,  then 

+00 
e2(t)    =     /JiKt)     "    *(t  -  x)|2dt 

+°°  +00 

=    /J*(t)|2dt   +    /J,Kt  -  T)|2dt 

+<x> 

-   /0o[^(t)^*(t  -  T)   +  l|>*(t)Tp(t  -  T)]dt 
+00 

=  2 /Ju(t)  |2dt  -  2Re[exp-j(jj  T[/u*(t)u(t  -  x)dt]] 

Since  the  first  part  of  the  result  is  proportional  to  the 
total  energy  of  the  signal  and  has  a  constant  value,  only  the 
second  term  is  going  to  make  £%t)    vary.   Therefore 


+00 


C(T)   =   /   U*(t)u(t  -  T)dt, 


*  =  complex 
conjugate 


is  defined  as  the  range  ambiguity  function;  only  c(t)  varies 
the  value  of  e2(t).   As  can  be  seen,  there  is  a  perfect 
identity  between  c(t)  and  the  autocorrelation  of  the  signal. 
The  ideal  case  would  be  to  have  the  ambiguity  function  with 
a  spike  at  the  origin  and  zero  anywhere  else.   That  would  mean 
that  the  only  situation  where  it  was  impossible  to  differen- 
tiate the  two  signals  was  when  they  were  on  top  of  each  other. 

As  an  example.  Fig.  4(a)  shows  the  autocorrelation 
function  of  a  pulse  of  duration  T,  Fig.  4(b)  its  uncertainty 
function  |c(x) |2.   In  part  (c)  there  is  a  rectangle  with  the 
same  height  of  c2(o)  and  with  the  same  area  under  |c(x)|2, 


24 
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-T 


T 


(b) 


— AT  - 


Fie.    4 


2  5 


that  is,  A  =  B.   As  it  can  be  seen,  the  base  of  that  rectangle 
is  a  measure  of  the  spread  of  the  curve  and  is  called  the 
delay  resolution  constant;  so 

_  /    lc(T)[2dT 

ZiT  —   -co — ' ' 

c2(o) 

sometimes  is  more  practical  to  define  At  as  a  function  of  the 

+00 

bandwidth.   Let   u(w)  =  F[u(t)]  =  fm   u(t)  exp  -  joit. 
since  by  Parsevall's  Theorem  F  [c(t)]  =  |u(w)|2 

2tt    /       |u(co)  I  Mw 


At    = 


[.OuHfdw]2 


if    the    effective    bandwidth    is    defined    as 

[/+c°|uU)|2dco] 

CO       1 


W 
e 


**iT\uM  Pdu 


At 


2  W 
e 


and  the  range  resolution  becomes 

AR  =  ^1- 
2 

All  these  definitions  make  sense  since  as  pointed  out 
before,  with  a  short  pulse  there  is  a  higher  ran-   resolution, 
but  to  have  a  short  pulse  there  is  need  of  a  wider  bandwidth, 
so  there  is  equivalence  in  the  two  expressions  for  At. 
3 .   The  Velocity  Ambiguity  Function 

Using  the  same  type  of  mathematical  logic  and  defining 

*(f)  =  FU(t)] 
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a  complex  correlation  function  is  determined  [5]  to  be 


K(fd)    =  /u*(2irf)u(2Tr£    -    2-nf£df 

=  /|u(t)|2    exp   j2fr£jtdt 

and  the  doppler  resolution  constant  [5] 

A£   =  /lK(fd)  l2dfd  =  /|u(t)  Tdt    E  J_ 
d       k2(o)       [/|u(t) |2dt]2"  Te 


where  Te  is  the  effective  duration. 

0 .      r  2vf o 

Since  ±a    =   

u     c 


df,  = 


2dv  fo 


which  infers   A f ^ 


2  f  o 


Av 


the  velocity  resolution  constant  is 
c 


Av  = 


2  fo  Te 


The  graphical  interpretation  for  Te  is  similar  to  the 
graphical  interpretation  given  for  At  in  Fig.  4  but  instead 
of  using  the  function  |c(t) |2,  |K(fj)  j2  is  used.   But  as  seen 
earlier,  the  optimization  of  one  resolution  implies  the  mini- 
mization of  the  other.   This  implies  the  need  to  study  a  two 
dimensional  correlation  function  in  order  to  see  the  mutual 
effects  between  the  two  resolution  parameters. 

So  if  the  transmitted  signal  is 

i|j(t)  =  u(t)  exp  j2iTfot 
the  received  signal  with  doppler  and  delay  will  be 
ip(t-x)  =  u(t-i)  exp  j27r(f0  -  fd)  (t-x) 
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The  mean  square  of  the  difference 

£  2  =  /  |iKt)  -  ij>  (t-x)  |  2dt 
will  yield  a  two  dimensional  correlation  function  of  the  form 

X(T,fd)  =  /u(t)  u*(t-T)  exp  j27Tfdtdt 

Easily  it  can  be  seen  that 

X(t,0)  =  c(t) 
X(0,fd)  =  k(fd) 

If  the  volume  under  |X(T,f,j)|2  is  determined  and 
divided  by  |X(0,0)|2,  an  equivalent  resolution  parameter 
called  the  effective  area  of  ambiguity  is  obtained. 

|X(0,0)|2 
It  is  easy  to  prove  [6]  that 

//|X  (T,fd)|2dTdfd  =  |X(0,0)|2 
A(i,fd)  =  1 

So  if  an  effort  is  made  to  improve  one  dimension,  the  other 
will  never  get  better.   But  given  one  acceptable  resolution 
in  one  dimension,  there  are  ways  to  find  the  optimal  solution 
for  the  other  dimension  if  X(x,fd)  is  known.   Fig.  5  is  an 
ambiguity  surface  for  a  monochromatic  pulse  with   a  smooth 
envelope.   It  would  be  a  reasonable  temptation  to  try  to 
design  the  waveform  from  a  desired  ambiguity  function.   That 
turns  out  to  be  impractical.   First,  a  criteria  to  design  the 
ambiguity  function  would  be  difficult  to  build.   Secondly,  and 
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most  importantly,  when  designing  a  radar  system  normally  there 
are  other  more  important  factors,  from  economic  to  space  and 
weight  constraints  as  well  as  technical,  that  dictate  the 
boundaries  for  the  waveform  that  we  have  to  pick.   Also,  there 
is  an  important  relation  between  the  chosen  waveform  and  the 
type  of  clutter  model,  and  the  approximation  of  the  clutter 
model.   Normally  some  flexibility  in  the  digital  processing 
is  desirable. 

C.   PULSE  COMPRESSION 
1 .   Concept 

Soon  it  was  found  that  trying  to  increase  the  maximum 
range  by  increasing  power,  keeping  the  same  pulse  width,  was 
the  difficult  way.   Isolation  problems  were  difficult  to 
handle  as  well  as  reliable  components.   The  only  way  to  in- 
crease the  energy  of  the  signal  was  by  increasing  the  pulse 
width.   But  an  increase  in  pulse  width  would  reduce  the  range 
resolution.   Since  the  real  problem  was  to  increase  the  time 
of  transmission  but  at  the  same  time  do  not  decrease  the 
bandwidth  (effective)  of  the  waveform  some  different  types  of 
waveforms  were  studied,  using  the  ambiguity  function  concept. 
The  processing  of  the  signal  should  be  such  that  the  trans- 
mitted waveform  with  high  time  of  transmission  and  small 
effective  bandwidth  (We)  should  be  converted  to  a  high  energy 
pulse  with  pulse  width  approximately  1/We.   This  is  called 
pulse  compression.   There  are  two  basic  ways  of  implementing 
pulse  compression  techniques.   The  active  way,  which  is  a 
time  domain  approach  which  basically  uses  correlators  in  the 
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detection  and  active  elements  to  produce  the  waveform  in  the 
transmission.   The  passive  way,  which  is  a  frequency  domain 
approach,  uses  passive  filters  to  generate  the  waveform  and 
matched  filters  in  the  receiver.   Each  of  the  two  has  some 
variations,  and  it  is  even  possible  to  build  a  system  as  a 
combination  of  the  two  processes.   Table  I  gives  a  relative 
performance  of  various  types  of  pulse  compression  techniques, 
2 .   Linear  FM  (Chirp)  Pulses  Compression 

First,  the  passive  generation  of  linear  FM  signals 
will  be  analyzed.   An  IF  pulse  generator  feeds  a  dispersive 
delay  filter  with  a  frequency  versus  time  characteristic  as 
in  Fig.  6(c).   The  signal  is  then  up  converted  and  only  one 
of  the  bands  is  transmitted.   At  the  receiver  the  signal  is 
down  converted  to  IF  and  then  passed  through  a  dispersive 
filter  with  the  opposite  slope  (Fig.  £    c) .    Using  at  recep- 
tion a  mixer  sideband  inverter,  the  same  type  of  dispersive 
filter  can  be  used  because  as  it  can  be  seen,  f~  and  f, 
(Fig.  6  b)   have  opposite  frequency  versus  time  character- 
istics.  In  fact,  with  this  technique  the  same  dispersive 
filter  for  both  transmission  and  reception  can  be  used. 

The  analysis  approach  to  this  type  of  processing  is 
the  following: 

Let  the  received  signal  be  of  the  form 

f(t)  =  exp  j  [  (a)0  +  u)d)  t  +  |ut2],  |t|  <  T/2 

=  0,  |t|  >  T/2 
where      u  =  ~  =    frequency  slope  of  the  dispersive  filter 
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(A)   GENERATOR  AND  DECODER 
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T  =  duration  of  the  transmit  pulse  envelope 
tod  =  doppler  shift. 

The  pulse  spectrum  is 

+  00 
F(u>)  =  f_m  f  (t)  exp  -jut  dt 

=  /«,  exp  j[(io0  +  wd  "  w)t  +  -^ y  t 2  ]  dt 

The  filter  transfer  function  is 

H(w)    =    exp    j  [(oo0    -    io)2/2u] 

which    implies    G(oo)    =    H(w)F(w)    =    exp    j  [  (io0    -    co)2/2]j] 

T/2  x 

/        exp    j  [  (o)0    +    03^    -    o))t    +    yyt2]    dt 
-T/2  * 

which    implies    the    time    function   at    the   output   of   the    filter 

is 

-1  1 

g(t)    =   F       [G(uO]    =    27  f  ~GM    exP    Jwtdw 

+00 
which    implies    g(t)    =   j^-  f_ m  exp   j  (to0    -    w)2/2y 

T/2  1 

/       exp   j  [  (io0    +    03d    "    W)T    +   tPt2] 
-T/2  Z 


dx    exp    jtot    dto 


Inverting    the   order    of    integration: 


T/2    +0°  ,  ,         fton  -  to)2 

(t)    =     /         /  exp    j  [(a)0  +  03d  -  u))t    +    2^T       +  2u  +   wt] 

-T/2   -°° 


dwdx 

T/2  1  u    2 

=     /     exp    j  [oo0t    +    wdT    +    7-yT2    +    y?-   ]    x 
-T/2  L  y 

+  OO  9  ,  . 

J"  exp    j  [-  cot  +    £-  +    tot    -    -^  o)j    do)      dx 
-co     ^     ^  L  2u  u 
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Call  v 


=  (ton  +  yr  -  yt) 
/T)T 


.2   _  WQ 


2y 


1       2  1       2 


-    yxt  -  w0t 


Multiplying  g(t)  by  exp  -v2  and  exp  v2 
1   TV  2  1 


+  0O 


Since , 


x    fm  exp  j  [to         /2y~  v]  2/2  y  doodx 

2/2~y    v  a) 


[a)  -  /2~y    v]  2    _    b^_ 
2y  2y 


+   v 


2y 


~—   +    v      -    —  a)    [  wn    +    yx    -    ytl 
2y  y         L     ° 


2  ,     w  . 

=     V        -     COT     +     o—    +     (Ot 


w0 


2y 


(0 


If  u  = 


CO 


2y    v 


/TU 


,  the  second  integral  will  be  of  the  form 


+  0O 


2y      /      exp    j  u     du    =    /2~y    vHT     exp    j  tt/4 


which    implies    g(t)    =     , exp   j  (co0t    -    j     yt 


1  ^2  TT-. 

+    4° 


4/2exP  jT(o)d    +    VJt)    dT 


Since    co0t    +    wdt    +    y    ^t2  "    y2    =    °°dT    +    ^Tt    "    2   yt2    +    (Jiot 


which    implies    g(t)    =     /— y-  cos  O0t    -    T  yt2)    sin    (a)d  +    yt)T/2 

/(^)  Z        (o>d  +  yt)/2 

/  9  T  2  i 

which  implies  g(t)  =  (-^  )cos  (<o0t  -  y  Vt2)    sin  x 


where 


x  =  (cod  +  yt)T/2 
x     2     2 
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sin  x 


From  these  results  it  can  be  seen  that  if  w,  =  0,  a 

a  '  x 

type   of   envelope    is    obtained  with   a   peak   at    the   origin   of 


magnitude    /r2AojT"T    .      So    the   higher   AtoT   the   better   the    output 

IT 

signal.   If  wj  /  0  there  is  a  shift  in  the  sin  x/x  curve  which 
implies  a  range  error  for  doppler  frequencies  greater  than 
zero . 

3 .   Matched-Filter  Approach 

The  matched-filter  approach  is  another  passive  method. 
The  analysis  of  the  former  will  give  us  a  basis  for  comparison 
Let  the  received  signal  be 

f(t)  =  cos  [(u>0t  +  wd)t  +  \   Pt2  ]  ,  |t|  <  T/2,0  elsewhere 

The  output  of  the  matched  filter  will  be  of  the  form 

g(x)  =  /*°°f(t)  h  (t  -  t)  dt 

where  h(t)   =Kf(-t)  and  F.[h(t)]  =  H(w)  is  the  filter  transfer 
function.   K  is  a  constant  factor  to  give  a  unity  gain.   In 
this  case,  K  =  /(-2yu  .   So 

TT 

g(T,03d)=  /^ZTk  ij^os  [(w0  +  wd)t  +  y  ^t2] 

TT 

cos     [oj0(t    -    t)     -    |y   (t    -    T)2]dt 
after   some   steps    [7] 


(T,<od)    -(Ijl  cos[(«n   ♦  ^)x]    sin    ("d*^      (T   -    |t|)] 


IT 

for       I t  !  <  T 


or,    g(T..d)    -   \    (T-   |t!)V^cos[(Wo  +  ^)t]    *i5_*    ,    |xj  <T 


where      x    =      u    2  —    (T    "     lTU 
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Since  x  is  not  a  linear  function  of  t,  the  shape  of 
the  curve,  unlike  in  the  previous  case,  is  not  a  sin  x/x  type 
of  curve.   Also,  since  w^x  does  not  represent  a  linear  shift 
in  the  axis  due  to  the  non-linearity  of  t  with  respect  to  x, 
if  wd  /  0  the  curves  become  distorted  and  not  symetric  with 
respect  to  the  vertical  axis  as  in  the  previous  case  (Fig.  7). 
The  advantage  is  of  course  the  maximization  of  the  signal  to 
noise  ratio. 

4 .   Use  of  Discrete  Frequency  Sequences  in  Pulse  Compression 

The  use  of  a  sequence  of  frequencies  each  with  a  pulse 
duration  of  t  in  order  to  get  a  good  autocorrelation  function 
is  a  very  common  method. 

Various  combinations,  from  linear  stepped  frequency 
to  randomly  chosen  frequencies,  are  used  in  order  to  increase 
the  main  lobe  with  respect  to  the  side  lobes  in  the  ambiguity 
function.   The  basic  format  is  that  of  Fig.  8 


A 


f, 


nr 


~t  • 


f. 


f 


f. 


Fig.  8 


where  the  waveform  of  the  nth  segment  can  be  expressed  as 


vn(t)    =   A   exp   j    2ir(£nt    +    <j>n) 
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The  mathematical  analysis  for  the  general  case  is 
quite  difficult  and  of  little  interest.   Since  only  a  few 
cases  are  of  interest,  consider  the  linear  stepped  frequency, 
that  is,  a  time  function  of  the  form, 


v(t) 


N-l 
n=0 


|n[u(t-ni)  -  u(t-(n  +  1)t)  ]  cos  (to   +  nAw)  t 


where 


Aoj 

N 


=  lowest  frequency  to  be  transmitted 

=  frequency  spacing 

=  number  of  frequencies  in  transmission, 


Assuming  that  there  is  coherency  between  all  frequen- 
cies, that  is,  cj)   =  0,  which  could  be  obtained  using  a  frequency 
synthesizer  with  a  master  oscillator,  the  matched  filter  would 
be  of  the  form  indicated  in  Fig.  9. 

nPut  I 


y 


r 


r 


t, 


A, 


Ah 


r  - 

r  ■ 

t 

If 

A 

a 

A 

VI 

it    summer 

output 


Fig.    9 
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If  xAf  >  1,  there  will  be  significant  range  ambigui- 
ties.  Also,  if  xAf  <  1  the  subpulse  filters  f   of  Pig.  9 

n      b 

will  overlap.   That  is  the  reason  why  in  most  cases  xAf  is 
made  equal  to  unity,  that  is,  x  =  ^.   With  that  assumption, 
the  convolution  of  the  impulse  response  of  the  matched  filter 
h(t)  =  v(-t)  will  yield  an  output 


e0(t) 


N-l 

nlQ   A*  exp  j  (a)o  +  nAw)t 


where    (N-l) t  <  t  <  Nt . 


If  A   =  1,  for  all  n,  then 
n  ' 


N-l 
eQ(t)  =  exp  j  ajQt  nS0(exp  j  Awt)n 


Since    | exp  j  Awt [  <  1,  then 


eft)  =  [exp  j  [w  +  (N-l)Ao)]  t]  x  sin  N  (Aoo/2)t 


o 


sin  (Aw/ 2) t 


Accounting  for  the  autocorrelation  of  the  rectangular 
pulse,  the  final  output  time  function  is 

e  ft)  =  (1-  |x|Af)  sin  N(Aw/2)t      r      (N-ljAio-, 
0  sin  (Ao)/2)t   C0S  LWo      2     J  L 

Fig.  10  shows  the  effects  of  xAf  f    1.   The  nules  of 
the  output  waveform  occur  at  sin  NAoo/2t  =  0. 
Which  implies  N  ~f-   t  =  ±  m  tt 

For  N  >  50,  the  first  side  lobe  is  13.46db  below  the  main 
lobe. 
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\      .  SUBPULSE AUTOCORRELATION 


Fig.    10 
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A  very  important  quantitative  factor  in  all  these 
waveforms  with  sharp  autocorrelation  function  is  the  pulse 
compression  ratio.   It  is  defined  as  the  ratio  of  peak  power 
after  compression  to  peak  power  before  compression.   Table  II 
gives  some  data  including  compression  ratios  for  some  typical 
linear  pre-pulse  generation  and  compression  circuits. 

Another  important  quantitative  factor  is  the  side  lobe 
reduction  factor  which  is  defined  as  the  power  ratio  of  the 
main  lobe  with  reference  to  any  side  lobe. 

Depending  on  the  specific  application,  it  is  sometimes 
preferable  to  sacrifice  the  main  lobe  peak  power  or  even  width 
(3  db  point) ,  but  get  a  higher  side  lobe  reduction.   In  the 
radar  problem  those  side  lobes  may  be  identified  with  targets, 
the  need  to  suppress  them  if  wanting  to  use  a  larger  dynamic 
range  is  obvious. 

There  are  various  types  of  weighting  functions.   The 
weighting  is  applied  in  the  amplitude  response  of  the  matched 
filter  in  order  to  deliberately  produce  a  mismatch. 

The  frequency  response  of  the  matched  filter  is  such 
that  the  signal  to  noise  ratio  is  maximized.   So  any  mismatch 
will  lower  that  signal  to  noise  ratio. 

The  loss  in  signal  to  noise  ratio  due  to  weighting 

is  called  the  loss  factor  L  .   It  is  defined  as 

rp  /xn    •  ,.  ,    r/rptoCt)  dt]2 
i   =  (S/N)  weighted  _  L  T  *■  J        i 

Ls    (S/N)  matched     Tf^z  (t)  dt 

w(t)  =  weighting  function 

T   =  processing  time  interval  or  the  length  of  the 
transmitted  waveform. 
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In  the  case  of  discrete  frequency  amplitude  weighting,  the 

loss  factor  is 

N 
[I    An]2 

Ls  =  -hr r8J 

N[E    A2! 

The  problem  of  weighting  is  in  effect  a  way  of  producing  a 
pulse  compressed  spectrum  that  yields  the  wanted  waveform. 

In  the  continuous  linear  FM  case,  the  spectrum  of  the 
match  filter  is  of  course  a  pulse  in  the  frequency  domain. 
Since  the  inverse  transform  of  a  pulse  is  a  sin  x/x  type  of 
function,  that  is  why  for  the  linear  FM  waveform  a  sin  x/x 
type  of  output  is  obtained.   From  the  Fourier  transform  pair 
it  is  known  that  a  wide  pulse  in  one  domain  corresponds  to  a 
sharp  pulse  in  the  opposite  domain.   So  the  amplitude  weight- 
ing functions  that  will  lower  the-  side  lobes  will  have  a  bell 
or  tapered  shape  characteristic  in  the  frequency  domain.   That 
is  why  the  time  output  of  a  matched  filter  with  a  pulse  in  the 
frequency  domain  as  transfer  function,  is 

1   w 
e„(t)  =  o—  /   a(<jo)  exp  j  go  t  d  co 
ov'    Ztt-w 

where  a(co)  is  the  weighting  function;  the  most  commonly  used 
weighting  functions  are: 
The  cosine  function  where 

a(uO  =  cos  j^     where  a(w)  =  1  at  the  center  of  the 

spectrum. 

So,  e0(t)    =   ^   /cos   gexpjutdw=    ^.^tw*        C0S    wt ' 
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The  main  lobe  will  be  between  nules.   The  nules  will 
occur  at  cos  wt  =  0 
which  implies   wt  =  (2n  +  1)  tt   ,  n  =  0,   ±1    ±2 

or    t  =  (2n  +  1)tt 
2w 

Considering  the  lobes  between  nules  the  side  lobe  reduction 
is  23. 5  db. 

Another  class  of  tapers  are  the  Hamming  functions 
which  have  the  general  form 

G  O)  =  a  +  (1  -  a)  cos  (  -^)   0  <  a  <  1 

where  the  center  of  the  spectrum  is  assumed  to  be  at  w  =  0. 

In  this  case  the  output  will  be 

^r-o_asinwt      1-a 
e°UJ  ~  2     wT~   2w(l-t*) 

for   a  =  . 54  the  side  lobe  reduction  is  42.8  db  which  represents 
a  much  larger  dynamic  range.   Table  III  gives  data  on  side 
lobe  suppression  on  frequency  coded  waveforms. 

D.   MOVING  TARGET  INDICATOR  (MTI) 

The  MTI  is  basically  a  processor  that  separates  moving 
targets  from  clutter.   Since  the  clutter  spectrum  is  a  very 
low  frequency  spectrum,  the  MTI  processor  must  be  a  high  pass 
filter.   In  the  early  MTI  processors  the  main  technical 
problems  were  in  the  processor  itself  clue  to  instability  of 
some  elements.   Today  the  main  problems  result  due  to  the 
existence  of  a  variety  of  types  of  clutter  which  make  the 
design  difficult.   If  only  the  elimination  of  one  type  of 
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clutter  is  required,  then  the  efficiency  of  the  MTT  is  good, 
but  if  various  types  of  clutter  must  be  handled  at  the  same 
time  by  the  same  processor,  the  efficiency  drops. 

Depending  on  the  parameter  used,  the  MTI  processors  are 
basically  divided  into  three  types:  those  that  use  the  phase 
information  of  the  returned  signal,  those  that  use  the  ampli- 
tude, and  those  that  use  both  phase  and  amplitude.   The  basic 
configuration  of  the  phase  processing  or  coherent  MTI  is  that 
of  Fig.  11.   The  use  of  two  oscillators,  the  stalo  and  coho , 
is  for  up  and  down  converting  the  waveforms  in  order  to  get 
a  perfect  phase  comparison,  since  at  low  frequencies  the  phase 
accuracy  measurement  increases.   Coherency  is  also  obtained. 

Consider  two  signal  returns  from  the  same  moving  target 
with  a  separation  in  time  by  T  =  1/ £R  where  f„  =  pulse 
rejection  frequency.   The  signals  presented  to  the  summer  are 

E-,  =  E  sin  (u)dt  +  <J)-,) 

E2    =    E    sin    [ojj(t  +  T)    +    4>2] 

*i  =  *z  =  4ttVa  =  * 

The   output    of   the  canceller   will    be 

Er   =    Ex    -    E2    =    2Esin(^)    cos[a>d(t    +   |)    +    <f>] 

R 

where  it  is  seen  that  the  cosine  waveform  is  modulated  by  a 

sin  ^S^-   .      So,  independently  of  t,  there  are  some  doppler 

iR 
frequencies  where  the  power  out  will  be  zero.   These  are  the 

frequencies  where  -^  =  niT .   That  is,   fd  =  n£R   (multiples 

of  the  pulse  repetition  frequency) . 
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■Ef  =  E(.)-E<t«-T) 


Fig.    11 
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Since  the  doppler  frequency  is  a  function  of  the  radial 

2Vr 
velocity  VR,  f,  =  — —  ,  there  are  also  blind  speeds  that 

correspond  to  the  blind  frequencies   f £ ,  =  nf„l 

1         v  d     R' 

Vblind  =  n(~2K   )  >  n  is  an  integer. 

Fig.  12  shows  the  power  response  for  a  single  delay  line 
MTI  processor.   As  it  is  seen,  the  multiples  of  the  PRF  are 
at  nules  but  the  other  frequencies  are  not  processed  in  the 
same  way.   That  is  very  far  away  from  the  ideal  case. 

The  clutter  spectrum  beating  with  the  pulse  repetition 
frequency  is  translated  to  every  multiple  of  the  PRF.   So  if 
the  clutter  spectrum  is,  after  translated,  of  the  form  indi- 
cated in  Fig.  13  (a)  the  ideal  MTI  processor  should  have  a 
response  as  in  Fig.  13  (b)  . 

As  it  is  known  the  single  delay  line  processor  response 
is  far  away  from  that  of  Fig.  13  (b)  .   One  way  to  improve  the 
MTI  processor  is  by  using  multiple  canceller  filters. 

For  an  n-stage  cascaded  canceller  (Fig.  14) 

(c~)n  =  2   sin   (— f  )    [9] 

b>l  XR 

where  S   and  S.  are  the  peak  output  power  and  peak  input  power 
The  problem  is  that  if  many  cascaded  blocks  are  used,  the 
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(sin  — jp4    function  may  cut  some  low  velocities  of  interest. 

The  transmitted  waveform  has  a  spectrum  as  in  Fig.  15. 
The  fact  that  there  is  more  than  one  line  in  the  spectrum  is 
of  no  importance  to  the  MTI  processor  since  all  the  lines  are 
at  blind  frequencies. 
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The  amplitude  processor,  also  called  noncoherent  processor, 
has  a  block  diagram  as  in  Fig.  16.   The  advantage  of  this 
processor  is  due  to  the  fact  that  stability  problems  almost 
disappear  since  there  is  no  coherency.   Since  most  high  power 
transmitters  are  not  very  stable,  this  is  a  real  advantage. 
The  big  disadvantage  is  that  it  needs  a  permanent  existence 
of  a  strong  clutter  return  since  without  it,  it  won't  work. 
Fig.  17  represents  the  phasor  diagram  that  explains  how  the 
noncoherent  processor  works.   Since  the  clutter  is  of  the 
form,  E   =  dc  value,  in  a  phasor  diagram  is  represented  by  a 
stationary  vector. 


Then 


Let    E    (t)    be    the    return   of   a  moving   target 
Eft)    =    IE    ft)  I    cos   co  At 


s  ^    J  '    s  ^    ;  '  d 

Es(t    +    T)    =    |Es(t)  |    cos     [aid(t    +    T)] 

Let  Wit    =    <j)    -    A(J>/2 

a)d(t    +    T)    =    $    +    A(J)/2    =    4>    -    A<J>/2    +    wdT 
which    implies    WjT    =    Ac}) 
Let    the    total    signal    at    time    t    (Fig.    17    b)    be 

^w  ■  tc  ♦  tsw 

and  the  total  signal  at  time  (t  +  T)  (Fig.  17  c)  be 
E^(t)  =  E1(t  +  T)  =  E^  +  E^Ct  +  T). 


From  Fig.  17  (d)  it  can  be  seen  that  due  to  E£  there  exists 
an  amplitude  difference  between  |E1(t)|  and  |E2(t)|.   This 
difference  in  amplitude  is  going  to  be  used  to  identify  the 
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E/t)=E(t+T) 


(b) 


(c) 


(d) 


EA=IEI-IEI 


moving  target.   In  order  to  get  better  results,  instead  of 
the  difference  in  amplitude,  the  difference  of  the  square  of 
the  amplitudes  will  be  used.   That  is  the  reason  for  a  square 
law  detector  in  Fig.  16. 
From  Fig.  17  (b) : 

Ef  =  E2  +  E|  -  2E2E|  cos  (<j>  -  Ac|>/2) 

(the  equation  relates  only  the  absolute  values  of 
the  vectors) 

From  Fig .  17  (c) : 

E22  =  E*  -  E|  -  2EcEs  cos  (4,  +  A(j>/2) 

(only  amplitude  relations) 

Then     E2  =  F2  -  E2  =  4E  E   sin  <f>  sin  ^ 
r     1     2      c  s      N      2 

since    (j>  =  ojjt  +  A  c|>  /  2 

and     Ad>  =  u  ,T  =  -^J^1  ,  T  =  ~ 
d      fR        £R 

which  implies 

E2  =  4EcEs  sin  ~^—   sin  (ojt  +  A<|)/2) 

So,  the  blind  speeds  are  the  same  as  in  the  coherent 

processor,  and  the  transfer  function  must  have  the  same  shape. 

But  if  E   =0,  that  implies  E   =0.   That  is  the  reason  why 

an  amplitude  processor  must  have  a  clutter  return,  and  the 

stronger  the  clutter  return  the  better  since  E   increases 
&  r 

linearly  with  E  . 
J  c 

Clutter  cancellation  can  also  be  made  at  IF  frequencies 
instead  of  at  the  video  part  of  the  receiver  (Fig.  18). 
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Let  the  two  IF  signals  separated  by  T  =  i-  be 

fR 

E1    =   E    sin    [2ir(flp   ±    f d)  t    -    <Pq] 

E2    =    E    sin    [27r(£IF    ±    £d)  (t    +    T)     -    <PQ) 

where        6      =    4tt£1FRo 
o        — 


Er    =   El  "  E2    =    2E    sin    [^  Cf  IF  -+  fd)T]    cos 

[2ir(£IF±£d)(t  +  ^)    -    *Q] 

is  the  output  of  the  summer. 

The  output  of  the  phase  detector  would  be 

EQ  =  E  sin  [Tr(fIF±fd)T]  cos  {2Tr[fdt  +  (fip  ±  fj)  |]  -  *q} 

where  the  video  envelope  is  of  the  form  E  sin  [it  (f  T„  ±  f ,)  Tl  . 

ir    d 

Since  in  an  MTI  processor  if  f ,  =  0  the  output  must  be  zero 
that  means  that  ^fjpT  =  niT  is  a  necessary  condition.   That  is, 
f   =  =-  =  nfR,  the  IF  frequency  must  be  a  multiple  of  the 
pulse  repetition  frequency. 

The  existence  of  blind  speeds  can  constitute  a  problem 
since  every  velocity  is  not  processed  in  the  same  way.   One 
of  the  most  often  used  methods  to  improve  that  situation  is 
with  the  use  of  staggered  PRF  systems.   There  are  some  quanti- 
tative factors  that  describe  the  efficiency  of  an  MTI  processor. 
Since  they  are  the  descriptors  of  the  MTI  processor,  it  is 
worthwhile  to  state  their  proper  definitions. 

Improvement  factor  is  defined  as  I  =  r0/r-;>  where  rQ  is 
the  output  target  to  clutter  ratio  and  r.  the  input  target 
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to  clutter  ratio.   This  definition  reflects  the  gain  as  well 
as  the  clutter  rejection  of  the  processor. 

Subclutter  visibility  is  defined  as  the  capacity  of  a 
radar  to  detect  moving  targets  in  a  clutter  environment.   A 
radar  with  x  db  of  SCV  is  one  that  is  able  to  detect  a  target 
over  a  clutter  that  has  a  signal  x  db  stronger  than  the  moving 
target.   The  SCV  cannot  be  used  as  a  parameter  for  comparison 
between  radars  since  the  target  to  clutter  ratio  is  a  function 
of  the  size  of  the  radar  resolution  cell. 

The  two  main  problems  with  MTI  processors  design  are: 
first,  a  correct  modeling  of  the  expected  type  of  clutter, 
and  second,  the  correct  shaping  of  the  transfer  function 
depending  on  the  type  of  clutter.   In  [2]  there  is  an  extensive 
description  of  the  different  types  of  clutter  and  modeling 
processes . 

When  the  platform  of  the  radar  is  moving,  the  stationary 
targets  will  have  a  doppler  return  different  from  zero.   This 
is  the  typical  case  of  the  Airborne  MTI  processors.   Various 
techniques  of  clutter  looking  and  automatic  tracking  must  be 
used.   The  techniques  get  much  more  complicated  and  the 
processor  becomes  much  more  dependent  for  specific  application. 

E.   DIGITAL  SIGNAL  PROCESSING,  DFT  AND  FFT 

Due  to  high  speed,  low  cost,  versatility  and  almost  inde- 
pendence of  external  condition  of  recent  digital  computers, 
the  processing  of  radar  signals  in  digital  form  covers  now  a 
wide  spectrum  of  the  radar  signal  processing. 
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The  two  basic  methods  of  approach  are  represented  in 
block  diagram  in  Fig.  19.   In  the  first  case  (Fig.  19  a)  the 
signal  is  first  sampled  and  quantized  in  such  a  way  that  the 
signal  is  converted  into  a  sequence  of  numbers.   The  numbers 
are  then  stored  and  the  processor,  using  arithmatic  and  logical 
operations,  manipulates  the  numbers  according  to  an  algorithm 
which  is  a  function  of  the  type  of  filter  we  want.   The 
resultant  numbers  are  then  dequantized  and  passed  through  a 
low  pass  filter  (D/A) .   The  second  method  (Fig.  19  b)  is 
basically  a  frequency  domain  approach  of  the  problem.   The 
signal  is  quantized  and  a  discrete  Fourier  transform  algorithm 
is  applied  to  the  quantized  signal.   The  spectrum  of  the  signal 
is  then  weighted  as  a  function  of  the  type  of  filtering  re- 
quired.  Then,  the  inverse  operations  take  place.   Due  to 
recent  improvements  in  faster  algorithms  to  implement  the  DFT 
this  second  method  is  becoming  widely  used.   The  theory  behind 
the  first  method  is  based  on  the  Z  transform  theory  [10].   The 
variable  Z  represents  a  delay  and  is  defined  through  its  Laplace 
transform  by 

Z  =  exp  ST 

where  1/T  is  the  sampling  frequency.   Given  an  analog  function 
f(t),  the  transform  of  f(t)  will  be 

F(Z)  =  nfQ  f(nT)Z"n 

The  two  basic  configurations  are  recursive  and  nonrecursive 
filters.   In  the  nonrecursive  case  the  output  is  of  the  form 

Eq  =  B^l  +  A1Z"1  +  .  ..  AnZ"n  ..]  , 
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where  the  coefficients  Ai  are  obtained  from  the  coefficients 
of  a  Fourier  series  expansion  of  the  frequency  domain  function 
of  the  desired  filter. 

The  feedback  or  recursive  filters  are  normally  obtained 
using  the  analog  transfer  function  of  the  required  filter  and 
substituting  S  by  a  chosen  relation  between  S  and  Z.   Usually 
the  bilinear  transformation 

[11]        S 


T        -1 

1  +  Z  L 

is  used  since  with  the  correct  sampling  frequency  it  is  very 
accurate  and  easy  to  implement.   Figure  20  shows  the  basic 
configurations  for  the  nonrecursive  (Fig.  20  a)  and  recursive 
(Fig.  20  b  )  filters. 

In  many  radar  applications  the  discrete  Fourier  transform 
is  being  used  especially  as  a  bank  of  filters  for  obtaining 
doppler  outputs.   The  understanding  of  what  it  is  and  its 
limitations  are  thus  fundamental. 

In  order  to  get  a  Fourier  transform  with  a  digital  pro- 
cessor not  only  the  input  waveform  must  be  in  digital  form 
but  also  the  values  of  the  Fourier  transform  must  be  repre- 
sented by  discrete  values.   Figure  21  gives  us  the  graphical 
derivation  of  the   discrete  Fourier  transform  pair.   The  input 
waveform  h(t)  (Fig.  21  a)  must  be  sampled  in  order  to  be 
manipulated  by  a  digital  processor.   The  sampling  operation 
is  the  same  as  multiplying  h(t)  by  a  string  of  delta  functions 
(Fig.  21  b) .   Since  multiplication  in  time  is  the  same  as 
convolution  in  the  frequency  domain,  the  resulting  spectrum 
is  that  of  Fig.  21  c)  .    Here  it  is  seen  that  if  the  samples 
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are  not  at  a  frequency  higher  than  the  Nyquist  rate  there 
will  be  an  aliasing  effect  that  will  distort  the  waveform. 
But  since  in  most  systems,  and  especially  in  any  real  time 
radar  system,  the  waveform  must  be  processed  in  a  finite  time, 
there  are  for  computation  only  a  finite  set  of  samples.   This 
is  equivalent  to  multiplying  h(t)A  (t)  (Fig.  21  c)  by  a 
rectangular  function  (Fig.  21  d)  .   Again,  the  previous  spec- 
trum must  be  convolved  with  the  spectrum  of  the  rectangle  to 
get  the  result  in  Fig.  21  (e) .   As  expected,  the  larger  the 

T  ,  the  smaller  is  the  distortion  in  the  spectrum.   But  as 
o '  L 

stated  before,  the  output  must  be  in  a  digital  format.   So  it 
is  necessary  to  multiply  in  the  frequency  domain  |H(f)  *  AQ(f) 
*  x(f)|,  by  a  string  of  delta  functions,  separated  by  1/TQ. 
That  is  equivalent  to  the  convolution  in  time  of  a  string  of 
impulses  separated  by  T   with  h(t)A  (t)  x(t).   The  final  result 
is  in  Fig.  17  (g) .   As  can  be  seen,  the  output  of  the  DFT 
processor  is  not  a  sampled  replica  of  |H(f)|.   Aliasing  effects 
and  finite  time  of  processing  are  responsible  for  the  amount 
of  distortion  introduced.   In  order  to  get  the  mathematical 
expression  for  the  DFT,  it  is  only  necessary  to  make  a  parallel 
derivation  to  that  of  the  graphical  derivation  already  dis- 
cussed.  The  time  domain  expression  for  the  output  (Fig.  17  g) 
will  be  obtained;  then  the  Fourier  coefficients  of  that 
periodic  function  represent  the  DFT. 

Let  h(t)  be  the  input  function.   The  operation  of  sampling 
h(t)  is  equivalent  to  a  multiplication  by 
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+  00 


A  ft)  =   E   6(t  -  KT) 
°       K=-oo 


T  00 

so   h(t)A  (t)  =   E   h  (KT)  6  (t  -  KT) 
°       K=-« 


Since    x(t) 


<  t  <  T   -  T/2 


0    otherwise 


is  to  be  multiplied  by  h(t)A  (t) 


N-l 


h(t)A  (t)x(t)  =   E  h(KT)  6  (t  -  KT) 
0  K=0 


where  NT  =  T 


o 


then  we  have  to  sample  the  Fourier  transform  of  h(t)A  (t)x(t) 
by  a  string  of  impulses  in  the  frequency  domain 


Mf>  =  ?~6(f  -  J/V 


-1  +oo 

which    infers      A^t)    =F     [A^f)]    =    TQ    TJ_m&(t    -    r   TQ) 


must   be    convolved  with   h(t)A    (t)x(t)    to    obtain   the    time   domain 
answer . 


hD(t) 


[h(t)A    ft)x(t)]    *    A-,(t)    = 


N-l  +  «> 

=    [    E    h(KT)6(t  -  KT)]    *     [T      rLra6(t-rT    )] 
K=0 

+oo        N-l 
=    Tn    r^  «>[    E   h(KT)6(t  -  KT-  rT    )] 
°         "       K=0 

Since  hn(t)  is  a  periodic  function  with  a  period  NT,  in 
order  to  avoid  time  aliasing  x(t)  is  chosen  in  such  a  way  that 
the  end  points  of  x(t)  don't  coincide  with  sampling  points. 
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If  they  did  coincide,  it  would  generate  an  additive  effect 
at  the  boundary.   But  NT  is  still  equal  to  T  .   Also,  the 
Fourier  transform  of  a  periodic  function  is  given  by  the 
Fourier  coefficients  of  a  series  expansion.   So 


Hn(nf   )    =   Hn(~)    =    E+°°  c      5(f-nf   ) 
D        oJ  D    i0  ~°°    n  o 


where 


,    T  -T/2 

1      ,   O    ' 


c   =  Fp—  /      hn(t)  exp  -  j  — f —  dt ,  n  =  0  ,   1 , .  . 
n     o  -T/2     U  Lo 


So ,  c 


n 


,    T  -T/2 

1  r  o     T 

T     .      0 

o  -T/2 


oo   n-1 


27mt 


£    ^   h(KT)6(t  -  KT  -  rT0)  exp  -jA^r 
r=-°°K=0  °        10 


dt 


Since  the  integration  is  only  over  one  period  TQ 


=  r 


T0-T/2 


n 


-T/2 


■N-1  ?   t ' 

E  S'(t  -  KT)  exp  -j—^ 
K>0  ° 


dt 


N-1 

Z  h(KT)  / 


T0-T/2 


n 


K=0 


T/2 


exp-j-2~   6(t  -  KT)  dt 


VurrT^       •  27rKnT 
E  h(KT)  exp  -j  — ~ — 

K=0  lo 


Since  T0  -  NT 


E  h(KT)  exp  -j  ^^  ,  n  =  0,    1,    2,. 


n 


K=0 


N 


and  the  Fourier  transform  of  hD(t)  is 


K—  u 


N 


NT- 
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n 


To  see  that  HD  (j^)  is  also  periodic,  it  is  necessary  to 
show  that  c^  is  periodic  with  a  period  NT. 

Let  n  =  r  +  N 


since  exp  -j 


2frK(r  +  N) 

N 


exp  -j 


27rKr 

N 


which  implies  c   (r  +  N)  =  c  (r) 
i       n       J  nK    J 


which  implies  HD  (^)  =  HD  (^) 


so  there  are  only  N  distinct  values  that  can  be  evaluated 
Normally  the  function 


n 


N-l 


H  (w)  =   Z   h  (KT)  exp  -j  2TrnK/N  ,  n  -  0,  1...  N-l 
1N  x  K=  0 

is  called  the  discrete  Fourier  transform  (DFT)  which  relates 
N  samples  in  the  time  domain  to  N  samples  in  the  frequency 
domain.   The  inverse  discrete  Fourier  transform,  is,  as  expected 

i  N_1  i      v 

h  (KT)  =  £  i    H  (Rnf)  exp  j  4jM  ,  K  =  0,  1...  N-l 

K-  0 
So  it  is  seen  that  with  very  simple  complex  multiplications 
and  additions  there  is  conversion  from  one  domain  to  the  other. 

The  Fast  Fourier  transform  is  an  algorithm  to  implement 
the  discrete  Fourier  transform  which,  by  reducing  the  number 
of  multiplications  reduces  the  processor  time  to  perform  the 
algorithm.   For  a  better  understanding  of  the  FFT,  it  must  be 
reduced  to  a  matrix  form. 


Consider  the  DFT 

N-l  9  v 

x(n)  =  T.      x  (K)  exp  -j       ,  n    0, 
K=0   ° 


N 


N-l 
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if  W  =    exp    -j    2tt/N,    then    in   a   compact    fo 
x(n)    =    WnK  x0(K) 


rra 


where      x (n) 


x(o) 


x(N-l) 


>    x0(K)    = 


x0(o) 


x0(N-l) 


W 


nK 


Wv 


Wv 


W 


(N-1)(N-1) 


the  examination  o£  the  matrix  equation  shows  that  there  will 

2 
be  N   possible  complex  multiplications  and  N(N-l)  complex 

additions.   If  N  =  2Y ,  and  the  FFT  algorithm  [12]  is  used, the 

number  of  multiplications  will  be  reduced  to  Ny/2  and  the 

number  of  additions  to  Ny .   Figure  22  gives  a  comparison 

between  the  FFT  algorithm  and  the  direct  computation  with 

respect  to  the  number  of  required  multiplications. 

As  it  was  shown,  complex  notation  is  greatly  used  in  the 
mathematical  structure  of  signal  processing.   So  the  need  to 
represent  signals  in  complex  notation  becomes  a  necessity  even 
in  their  implementation.   It  will  be  seen  how  a  sine  wave  type 
of  function  is  converted  in  order  to  be  represented  by  a 
complex  number. 

In  Fig.  23,  let  the  input  signal  be 


S(t)  =    A   sin    [0Ip   ±    ud)t    +    <fr0] 


where  wTT.  and  go,  can  be  thought  as  the  intermediate  frequency 
IF       d 

and  doppler  frequency  in  a  radar  receiver.   If  the  local 
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oscillator  (I/Q)  generates  2  cos  ooTp  and.  -2  sin  to   t  which 
can  be  obtained  with  a  90°  phase  shifter,  the  in  phase  channel 
(I)  and  the  out  of  phase  channel  (Q)  outputs  are  of  the  form 

I  =  A  cos  (to  -jt  +  cf)  ) 

Q  =  A  sin  Odt  +  <|>o) 

Since  the  low  pass  filter  retains  only  the  difference  frequency 
The  I  and  Q  signals  can  then  be  thought  as  the  real  and  imagin- 
ary part  of  the  complex  signal 

Z  =  I  +  j  Q 

If  the  I  and  the  Q  channels  are  sampled  and  converted  to 
sequences  of  numbers  x„  and  y„,  the  resultant  sampled  pair 
can  be  considered  a  complex  digital  word. 

AK  =  XK  +  ^  yK  =  lAK'eXp  ^  *K 
where   x„  =  A  cos  (to  ,  K  T  +  <j>  ) 

y   =  A  cos  (to  ,  KT  +  A) 

'  n  d      H  o 

|AK|  =  A,  cDK  =  todKT  +  ^ 

assuming  a  sampling  frequency  f  =  ^.  Since  t^  varies  linearly 

with  K,  the  complex  digitized  signal  can  be  interjected  as  a 

vector  of  amplitude  A  that  shifts  at  every  sampling  time 
KT,  by  todKAT. 
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F.   SYNTHETIC  APERTURE  RADAR 

Using  the  pulse  compression  techniques  already  mentioned, 
it  is  possible  to  achieve  very  small  values  of  range  resolu- 
tion.  In  many  radar  applications  such  as  aerial  photography 
not  only  a  very  small  range  resolution  is  necessary  but  also 
an  equivalent  small  aximuth  resulution  is  desired.   The  con- 
ventional method  of  reducing  azimuth  resolution  by  very  small 
beamwidths  will  always  require  big  antennas  which,  especially 
in  airborne  radars,  is  a  problem  since  space  is  a  limiting 
factor.   Also,  the  fact  that  azimuth  resulution  is  an  increas- 
ing function  with  distance  makes  it  many  orders  of  magnitude 
higher  than  the  range  resolution  even  for  small  distances. 

Consider  a  radar  without  a  pulse  compression  processor. 
Then  the  range  resolution  6R  is 

6R  =  j    p  c  =  speed  of  propagation 

B   =  bandwidth 
w 

If  the  antenna  has  an  aperture  D,  the  beamwidth  is  6g  =  -p-  , 

So  the  azimuth  resolution  6.-  is 

x      Dn   _  RA      R  =  distance  to  the  target 

6AZ  =  RGB  "  T 

if  B   =10  MHZ,  f  =  5000  MHZ,   D  =  10  yds  and  R  =  100  mi, 

w 

then     SR  =  16.4  yds 
6    =  1300  yds 

Synthetic  aperture  radar  is  a  way  of  processing  radar 
signals  in  order  to  get  an  equivalent  antenna  aperture  much 
bigger  than  the  real  antenna  to  reduce  6AZ>   In  the  synthetic 
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aperture  case  a  single  antenna  is  translate   i 

^  Liansiated  along  a  line 

The  received  signals  are  then  stored,  and  processed  afte/ 
the  radiating  element  has  travelled  &   ^   ^^   ^ 

synthetic  aperture  concent  ran  ho   • 

oncept  can  be  viewed  through  two  different 

aspects.   From  a  doppler  viewpoint  or  from  »  t 

F      or  trom  a  linear  array  view- 
point, the  basic  principle  nf  tv,„  a 

r-ixiiLipxe  or  tne  dopnler  PffPrt  „• 

^-ttI  crrect  viewpoint  is 

that  there  exists  a  one  to  one  eorrespondence  between  the 
!ong  track  coordinates  of  a  reacting  object  and  the  instan- 
taneous  doppler  shift. 

Consider  Fig.  24  in  wMch  a„  airplane  ^^  ^  ^  ^.^ 

h0  wxth  a  constant  speed  v  and  with  a  radar  antenna  that 

illuminates  the  area  A.   The  ra^r  nn  +  Q 

me  radar  antenna  does  not  rotate 

but  simply  moves  with  the  airplane. 

If  there  is  no  pulse  compression  and  the  width  of  the 

transmitted  pulse  iq  T   i  +■  a  c    i 

PUise  is  T,  lt  is  known  that  the  slant  resolution 

constant  is 
6 


'R   2  s  2~T 


w 


so  the  corresponding  ground  resolution  constant  is 
PR  =  SR   sec  if,  =  £L  sec  <(, 

Also,  from  Fig.  24  it  is  seen  that  L  .  BR  where  B  is  the  b 
width  of  the  transmitting  antenna  and  Le£f  is  the  the  along 

track  resolution  p  . 

x 


earn 


Since 


D 


then 


px  =  Leff  =  "d-  R 
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and  as  pointed  out,  it  is  impossible  to  increase  D  beyond 
certain  limits  in  order  to  decrease  p  .   But  if  a  relation 
is  found  between  doppler  frequency  and  the  position  x  of  the 
radar,  then  it  is  possible  to  make  an  azimuth  discrimination. 

From  Fig.  25  it  is  seen  that  the  radial  velocity  with 
respect  to  the  target  is 

v   =  v  sin  6 

So  the  doppler  shift  will  be  f ,  =  ~  f   sin  0  where  f   is  the 

a         c   o  o 

transmitted  frequency.   But  if  6  is  small,  which  is  the  normal 
physical  situation  since  the  total  illuminated  area  is  small 
and  R  is  big,  then  sin  8  -  6  =  x   -x,  which  makes 

±  ,   -   T-=r   (X     "   Xj 

d    AR    o     ' 

So  it  is  seen  that  with  a  frequency  analysis  of  the  received 
signal,  it  is  possible  to  get  an  azimuth  resolution  that  is 
only  a  function  of  discriminating  between  the  closeness  of 
two  frequencies.   The  basic  principle  for  the  explanation  of 
synthetic  aperture  radar  using  the  linear  array  theory  is  the 
following : 

In  a  linear  array  various  elements  are  fed  at  the  same 
time  and  the  received  echo  is  also  received  at  the  same  time, 
but  due  to  path  differences  a  receiving  pattern  is  generated. 
In  the  synthetic  aperture  case  there  is  only  one  transmitting 
element  that  moves  with  constant  speed  along  a  line.   So  a 
long  antenna  will  be  formed  not  by  physical  means  but  by  signal 
processing.   After  the  radiating  element  has  travelled  a  dis- 
tance L  r  r-   corresponding  to  the  illumination  time  of  the  real 
ef  f       1  b 
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beam,  the  stored  signals  after  weighting  and  phase  shifting 
resemble  the  signals  received  by  a  big  linear  array. 

Two  cases  are  important  to  differentiate:  the  unfocused 
synthetic  aperture  and  the  focused  synthetic  aperture. 

In  the  focused  synthetic  aperture  radar  if  e  =  0  (Fig.  26  a) 
by  signal  processing  (for  R  f   «)  or  by  real  situation  (R  =  °°) 
the  signals  all  arrive  with  the  same  phase.   In  a  real  linear 
array  the  angle  selectivity  is  provided  only  during  the  recep- 
tion of  the  signals.   But  in  the  synthetic  case,  since  only 
one  element  is  radiating  and  is  moving,  it  is  necessary  to 
account  for  the  phase  shift  due  to  the  transmitted  and  received 
path.   So,  due  to  the  fact  that  all  the  round-trip  phase  counts 
for  the  formation  of  the  receiving  pattern,  as  opposed  to  the 
real  case  in  which  only  the  received  path  is  responsible  for 
the  receiving  pattern,  the  effective  beamwidth  is 

3  x.c  -  T~i instead  of 


" 


ef£    2  Leff  Leff 

Since  the  length  of  the  synthetic  aperture  radar  is  equal 
to  the  distance  corresponding  to  the  time  of  illumination, 
L    =  ^r  ,   where  D  is  the  horizontal  aperture  of  the  physical 
antenna  and  R  the  distance  to  the  target.  (Fig.  24),  then  the 
azimuth  resolution  is 

6AZ  =  3eff  R  =  2  Le££  R  =  2 

The  azimuth  resolution  is  independent  of  X    and  R  in  the 
focused  case  which  is  of  course  an  important  result.   This 
means  that  the  smaller  the  antenna  the  higher  the  azimuth 

resolution. 
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In  the  unfocused  case  it  is  important  to  account  for  the 
phase  difference  between  the  center  and  the  ends  of  the  equiv- 
alent linear  array.   A  quantitative  interpretation  of  the 
unfocused  synthetic  aperture  radar  can  be  explained  as  follows 

Consider  Fig.  26  (a):   Due  to  the  fact  that  R   /  R,  there 

o     ' 

is  a  phase  shift  between  the  center  and  the  end.   The  total 
voltage  is  the  vector  sum  of  small  discrete  voltages  that  are 
received  along  the  path  L  ff  (Fig.  26  b) .   If  the  number  of 
discrete  voltages  is  large,  it  is  possible  to  make  the  approx- 
imation of  Fig.  26  (c)  and  determine  the  relationship  between 
the  resultant  voltage  V  and  the  linear  sum  of  the  discrete 
voltages  V  .   If  the  3  db  point  is  defined  as  the  boundary 
for  the  drop  in  voltage,  a  must  be  approximately  tt/2.   But, 
when  a  -  tt/2,  the  radius  r,  (Fig.  26  c)  is  perpendicular  to 
the  tangent  t-.  ,  implying  that  the  phase  difference  between 
the  voltage  at  the  end  and  at  the  center  must  approximate  45°. 

This  sets  a  limit  on  L  rr.   So  when  R  =  R   +  A/8  (round  trip 

eff  o       v 

<  >  A/4)  the  conditions  are  met.   Using  Fig.  26  (a): 
(Lef£/2)2  +  R2  =  (R  +  A/8)2  =  R2  +  ^  +  |i 

\   2  

if  ^-j  is  very  small  compared  to  the  other  terms  L  ££  =  /  RA 


Since  07J1  =  .455  , m   -~    /  A~ 

3db        Le£f    2    -T 


the  azimuth  resolution  would  be  6^    -   R0  -  j   /XR 

Since  6A7  is  a  function  of  R  '  ,  this  is  a  solution  between 
the  focused  case  and  the  conventional  case.   Figure  27  gives 
a  relation  between  the  three  cases. 
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The  signal  processing  theory  for  synthetic  aperture  radar 
can  be  understood  through  the  analysis  of  the  ambiguity  func- 
tion of  the  received  signal.   If  it  is  possible  to  separate 
the  range  and  azimuth  as  the  two  ambiguity  factors  from  the 
general  ambiguity  function,  then  the  determination  of  the 
azimuth  ambiguity  function  and  consequently  the  azimuth  reso- 
lution function  are  immediate. 

Consider  a  transmitted  signal  f(t).   After  reflection  in 
the  ground  or  target  with  reflectivity  p(x,y,z),  the  received 
signal  will  be 

S(t)  =  ///  p(x,y,z)  f  (t  -  ^)  dxdydt  [13] 

where  the  integration  is  made  over  the  illuminated  region, 
and  R  is  the  distance  between  (x,y,z)  and  the  radar  antenna. 
The  processing  of  the  signals  in  fact  recovers  a  signal  pro- 
portional to  p(x,y,z).  Using  matched  filter  techniques,  the 
output  will  be  of  the  form, 

eQ(t)  =  /  f*(t  -  -~4s(t)  dt 

=  ////  P(x,y,z)  f  (t-^)  f*  (t  -^-)  dtdxdydz 

where  R'  is  the  distance  from  the  antenna  to  the  point 
(x',y',z')  corresponding  to  ptx'jy'jZ1).   Let's  find  the 
ambiguity  function  of  f(t).   By  definition  it  will  be 

<Kx,y,z,x',y',z')  =  /  £(t  -  ^)  f*  (t  -  2|1)  dt 

where  *  means  complex  conjugate. 
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So 

eQ(t)  =  fffftyp   (x,y,z)  dxdydz, 

is  a  weighted  average  of  p(x,y,z)  where  the  weighting  function 
is  i^(x,y,z,x'  ,y'  ,z')  . 

Let     f(t)  =  g(t)  exp  j  uiot 

then        i    r       r+         2R^   *  r  +         2R'.      .    ,2R    2R\  ,+ 
*  =  f   g(t  -  — )  g*  (t  -  — )  exp  3  w0(—  -  — )  dt 

Assume  the  energy  is  transmitted  in  finite  amounts  of 
time.   If  during  those  finite  lengths  of  time  the  exponential 
does  not  vary  too  much,  then  during  a  transmission  the  expon- 
ential part  can  be  considered  constant  although  varying  between 
transmissions.   So  it  can  be  taken  out  of  the  integral.   The 
result  will  be 

v  r  r         f*  2^   *  r+  2R'    ,.  ,        .     r2R     2R' 

ij>  =  Z  [  /  g:(t  -  — )  g*  (t — )  dt]  exp  -j  u)o(- — ) 

If  gg*  is  constant  during  a  transmission  period,  it  can  be 
moved  out  of  the  summation  and 

i>    =    [/gg*dt]  2  exp  -j^0[-—  -  — ] 

It  is  obvious  that  the  first  factor  is  the  range  ambiguity 

function,  and  the  second,  the  azimuth  ambiguity  function. 

So,  the  azimuth  ambiguity  function  is 

r2R     2R' 
ij,AZ  =  Z  exp  -3    u)Q(- — ) 

Figure  28  indicates  a  real  situation.   From  that  figure 
it  can  be  seen  that: 
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Fig.     2, 


X  =  vt 


R  =  /R  <  + 


17 


Ro  +  2TT 
o 


Ri  =  /R-  +  (X  -  X') 


I  ^  2 


R   +   U  -  X') 

o      2R 
o 


if  x  <<  Rq  and  (x2  -  x')  <<  Rq ,  which  is  basically  the  same 
restriction  given  to  6  in  Fig.  25. 

If  the  radar  transmits  with  a  PRR  of  fD  =  1/T,  then  x  = 

K 

nvT   are    the   positions   where    transmission   occurs. 


which    infers      \\> .      =    z    exp    -j    w    (—  -    -^j 

AZ  r       J       o v  c  c 


=    Z    exp 


2^   (2xx-    -(x')2) 
o 


N/2 


=    exp    +    j    ^ — J—     z      exp    -j    4tt(x'AR   )  nvT 

-N/2  ° 


c  R. 


where   L,    the    synthetic   aperture   length    is 

L   =   NvT    ,        v   =    aircraft    speed 


which    infers    t|/.7 


pXn    i    ^o(x')2      sin    (N  +  1)4ttx'vT/2ARq 

p    J      Re  sin    4ttx*vT/2AR 

o  '         o 


f  x  '  ")  2 

where  exp  j  ~-^ is  only  a  phase  factor.   The 

Roc 
sin  (N +  l)a/sin  a  is  the  amplitude  factor  of  interest.   The 

azimuth  resolution  is  given  by  the  3  db  points  separation. 


If 


N 


N 


^-r-  -    1,    the    3  db     point    is    at 


2ttx'  L    _    ,     . 
AR„ 


and 


s       =   2x,   =  i^ARa 

AZ  ttL 
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Since   L  =  3RQ  =  ^  RQ 

6az  =  ~t~  D  K  2 

which  is  a  result  in  agreement  with  that  already  derived. 

In  generating  the  synthetic  aperture  radar,  in  order  to 
reduce  the  side  lobes  of  ^A7,  it  is  sometimes  useful  to  weight 
the  returned  signals  before  combining  them.   In  the  case  of 
the  focused  antenna,  phase  compensation  must  also  be  used. 
So,  if  N  returns  are  processed  in  order  to  get 

L  =  N  v  T 

the  focused  processing  will  be  of  the  form 

I  S  [exp  j  \b    ]W 
nL   ^  J    ynJ  n 

and  the  unfocused  processing  of  the  form 

Z  S   W 
n   n 

where  W   is  the  weighting  function  and  S^  are  the  discrete 
n  n 

returns.   This  is  the  basic  processing  technique.   Various 
methods,  from  optical,  electronic,  as  well  as  acoustic,  are 
used  to  generate  the  synthetic  aperture  radar. 


V.   RECENT  DEVELOPMENTS  IN  TWO  MAJOR  AREAS 

A.   DIGITAL  MTI 

The  structure  of  an  MTI  processor  clearly  points  into  the 
direction  of  a  digital  implementation.   The  increase  in  digital 
data  rates  associated  with  the  cost  decay  of  present  digital 
components  led  to  the  implementation  of  reliable,  flexible, 
low  cost  digital  MTI  processors. 

Fig.  29  is  a  typical  block  diagram  of  one  of  the  types  of 
MTI  digital  implementation.   I  and  Q  channels  are  used  in 
order  not  to  lose  3  db  on  the  average  due  to  blind  phases. 
The  structure  is  parallel  to  the  analog  MTI  processor.   The 
use  of  A/D's  and  D/A's,  storage  devices  instead  of  delay  lines, 
make  up  the  differences  with  respect  to  the  analog  processor. 
The  capacity  of  the  memory  depends  basically  on  the  number  of 
range  cells  and  on  the  use  of  multiple  cancellers. 

The  quantization  process  introduces  a  new  dimension  to  the 
noise  problem:  that  is,  the  quantization  noise.   Quantization 
noise  is  present  along  all  the  dynamic  range  of  the  processor, 
but  the  errors  at  the  extremes  are  of  particular  importance 
to  the  MTI  processor.   At  the  lower  level,  clutter  cannot  be 
cancelled  below  one  bit  of  quantization;  at  the  upper  level, 
abrupt  clipping  distorts  the  signal  creating  additional  noise. 
This  extra  noise  affects  the  improvement  factor  (I)  of  the 
processor.   If  it  is  assumed  that  signal  and  quantization 
noise  are  uncorrelated ,  then  for  a  single  canceller,  the 
improvement  factor  is  [14], 
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Il 


1    +    (opVE2)  [P,(T) 


_es. 


(T) 


l  -  p  (T)  .+   Ca77Pi^Tl~T~pTn7T   ' 


where 


p£  =  correlation  coefficient  of  quantization  error  for 
clutter  alone 

pes=  correlation  coefficient  of  quantization  error  when 


both  clutter  and  signal  are  present 
a   2  -   noise  power 
E   =  signal  voltage 
P.   =  input  power  clutter 
p(T)  =  correlation  coefficient  of  clutter  signals. 


Since  o  2  /E2  =  0   and   p  (T)  =  0 
e  e 


1  -  P(T)  +  a^/Pi( 


which  differs  from  the  equivalent  analog  factor  only  by  the 
term  0   2/P-j^. 

Also,  for  a  double  canceller,  the  improvement  factor  I 
is   [14] 

1 


I 


1  "  j   P(T)  +  ~   p(2T)  +  a£2/Pic 


For  the  triple  and  higher  order  cancellers,  the  expressions 
get  more  complicated  but  the  derivation  is  the  same  for  all 
cases.   Also  in  I  ,  the  only  difference  from  the  analog 
expression  is  the  factor  a2/P-  .   Assuming  a  uniform  error 

r  E    IC  ■ 


distribution 


a  2  = 


F  2 


,n- 


G    12(2il_-  1/2)  2 
..here  Em  =  saturation  voltage  of  the  digital  register 
n  =  number  of  bits  used  in  the  quantization. 


Pic/ae2  is  called  the  improvement  factor  limitation. 

Figure  30  relates  the  improvement  factor  limitation  with  the 

total  number  of  bits  for  different  values  of  P.  /E  2  . 

ic   m 

In  many  applications  the  fact  that  the  MTI  processor  has 
blind  speeds  at  multiples  of  the  PRF  is  completely  undesirable. 
The  method  usually  used  to  attenuate  this  effect  makes  use  of 
staggered  PRF.  The  scan- to- scan  staggered  PRF  is  less  efficient 
but  requires  less  complex  hardware  than  the  pulse- to- pulse 
stagger,  but  for  some  applications  the  first  one  is  inadequate. 
The  staggered  PRF  technique  is  based  on  the  use  of  different 
spacing  between  transmitted  pulses  which,  at  the  receiver,  are 
properly  delayed  (de- staggered) ,  such  that  the  spacings  at 
the  input  of  the  delay  line  cancellers  are  all  equal.   The 
necessary  use  of  delays  makes  staggered  PRF  perfectly  matched 
for  digital  techniques.   Since  the  concept  is  completely 
defined,  the  recent  trend  has  been  focused  in  the  development 
of  optimization  techniques  in  order  to  improve  the  performance 
of  the  processor.   The  first  difficulty  with  any  optimization 
problem  is,  of  course,  the  definition  of  optimum.   Some  try 
to  optimize  the  improvement  factor,  others  to  optimize  the 
signal  to  clutter  gain  (SCG)  within  a  finite  number  of  fre- 
quency slots,  or  even  the  optimization  of  some  defined  indicator 

An  important  factor  in  the  optimization  process  is  the 
statistical  model  that  is  chosen  for  the  clutter  power  dis- 
tribution.  It  has  been  verified  that  a  Gaussian  distribution 
centered  at  zero  doppler  frequency  is  adequate  in  most  cases 
as  a  power  density  function  of  clutter.   Various  standard 
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deviations  have  been  calculated  for  different  types  of 
clutter  [2] . 

Figure  31  is  a  block  diagram  of  a  staggered  PRF  processor. 
The  pulse  staggered  sequence  is  such  that  the  i^'1  pulse  is 
delayed  (i  -  1)T  +  AT.  seconds  with  respect  to  the  first  pulse. 
In  the  processor,  after  de- staggering ,  the  pulses  are  weighted 
(W.'s  in  Fig.  31).   By  observation,  the  impulse  response  of 
the  processor  is 

h(t)  =  W1  6(t  -  AT.,)  +  W2  6[t  -  (T  +  AT2)] 

+  ...  +  WN  6{t  -  [(N-l)T  +  ATN]} 

which  implies  a  transfer  function 

H(u>)  =  il1    \ij    exp  {-ju  [(i  -  1)T  +  AT.]} 

and  a  filter  power  response 

G(co)  =  H(oo)  H*(uO  ,      *  =  complex  conjugate. 

With  G(w)  and  the  power  spectral  density  of  the  clutter, 

different  optimization  criteria  can  now  be  devised.   One  factor 

is  common  to  all  techniques.   The  A-'s,  W.'s  and  the  clutter 

n  i   '   i 

variance,  assuming  a  Gaussian  distribution,  are  the  parameters 
that  will  optimize  the  processor, given  an  optimization  pro- 
cedure.  In  the  analysis,  it  is  always  assumed  a  step  scanned 
array  antenna  in  order  that  for  each  position,  N  pulses  are 
transmitted.   Figure  31  is  a  processor  for  a  given  range  cell. 
If  there  are  K  large  cells,  the  general  block  diagram  is 
presented  on  Fig.  32,  where  W.'s  and  the  summing  network  is 

presented  for  the  K.  cell. 
^  i 
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The  ideal  MTI  filter  is  a  high-pass  filter  with  rejection 
from  zero  frequency  to  the  maximum  frequency  of  the  clutter. 
Rejection  may  also  occur  after  the  maximum  desired  doppler 
frequency  but  it  is  not  absolutely  necessary. 

For  a  better  understanding  of  staggered  PRF  some  general 
concepts  should  be  explored.   If  a  constant  PRF  is  used,  not 
only  the  dc  component  of  the  doppler  spectrum  but  the  whole 
clutter  spectrum  is  translated  to  multiples  of  the  PRF.   It 
is  shown  [15]  that  with  staggered  PRF  the  clutter  spectrum 
still  is  translated  to  multiples  of  the   average  PRF,  but  the 
spectrum  of  the  signal  is  dispersed  into  N  (number  of  A-'s) 
separate  frequency  lines.   So,  if  a  target  produces  a  doppler 
frequency  that  is  a  multiple  of  the  basic  PRF,  only  one  part 
in  the  whole  signal  is  distorted  by  the  clutter  spectrum 
where  the  remaining  M-l  parts  of  the  signal  may  possess  enough 
power  to  be  detected.   An  MTI  filter  which  notches  in  a  region 
around  dc  as  well  as  in  multiples  of  PRF  can  eliminate  the 
clutter,  but  if  staggered  PRF  is  used  only  one  part  in  N  of 
the  signal  is  lost,  thus  increasing  the  probability  of  detec- 
tion of  targets  with  potential  blind  velocities. 

In  [16]  an  optimization  process,  based  on  a  performance 

index  P  defined  as 

Bu     1 
P  =  /» u   df 

Bl  [S*(f)]M 

is  developed.   B-.  and  B   are  the  lower  and  upper  bounds  of  the 
velocity  region,  S2(f)  is  the  signal  to  clutter  gain  and  M  is 
a  parameter  that  reflects  the  emphasis  the  criteria  puts  on 

92 


spectral  amplitude  variations.   As  M  increases,  smaller 
spectral  amplitudes  have  more  influence  in  the  final  result 
of  the  integral. 

For  an  understanding  of  the  significance  of  P,  it  should 
be  noted  that  good  clutter  attenuation  as  well  as  a  high  im- 
provement factor  correspond  to  a  good  performance  index. 
Using  the  optimization  procedure  described  in  [16]  and  using 
M  =  8,  practical  results  (figures  32,  33,  34)  for  different 
values  of  N  and  (B   -  B,)  were  obtained.   It  can  be  seen  that 
as  N  increases,  the  minimum  value  for  the  signal-to-clutter 
gain  in  the  velocity  region  increases  as  well  as  the  improve- 


ment factor  I  =  SCR(f)  (Fig.  32,  33),  but  with  a  fixed  N  and 
the  velocity  region  increased  on  the  upper  bound,  the  minimum 
value  of  SCR(f)  in  the  velocity  region  decreases  (Fig.  33,  34) 

As  was  said  earlier,  the  clutter  variance  a  is  also  an 
important  parameter  in  the  optimization  process.   Fig.  35 
relates  a  with  the  improvement  factor  for  the  case  of  N  =  5 
and  B   =  10  B, .   In  the  solid  line  the  interpulse  periods 
T  +  A-  were  optimized  for  the  value  of  o  =  0.02  B,  and  the 
coefficients  W.  were  optimized  in  parallel  with  the  variation 
of  a  in  the  horizontal  axis.   On  the  dotted  line,  both  coef- 
ficients and  interpulses  were  optimized  for  a  =  0.02  B-,  ,  and 
the  variations  of  I  with  a  were  recorded.   This  situation 
occurs  when  the  real  clutter  parameters  are  not  those  used 
in  the  design  of  the  filter.   It  is  apparent  that  there  are 
not  too  many  differences  between  the  two  cases. 
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SCR  gain  as  a  function  of  Doppler  frequency; 
N  —  S.B  ~  4 Off,  .  a  -  0.028 ,  .Optimum  intervals:  1.007. 
1.021,  1.000,  1.209.  Optimum  coefficients:  1.000, 
—3.784.  5.432.  -3.273,  0.626. 
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Improvement  factor  as  a  function  of  clutter  spectral  width 
o;N  ~B,  Bu=  1 0ff,  .  Coefficients  for  dotted  curve  (optimum  for  a  = 
0.02£,):  1.000,  -4.258,  6.242,  -4.031,  1.048.  Intervals  for  both 
tolid  end  dotted  curves  (optimum  for  o  =  0.025,):  1.111,  1.000, 
1.091,1.058  (c.f..  Fig.  8). 
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Another  criterion  to  optimize  the  staggered  PRF  processor 
is  developed  in  [17].   Instead  of  maximizing  the  average  sig- 
nal-to-clutter gain  over  the  entire  velocity  region,  the 
signal-to-clutter  gain  is  maximized  in  small  Af  intervals 
within  the  velocity  region.   The  number  of  intervals  is,  of 
course,  directly  related  to  the  complexity  of  the  processor. 
Using  the  optimization  algorithm  developed  in  [19]  ,  a  final 
block  diagram  for  the  processor  is  reached  (Fig.  36,  37).   In 
Fig.  36  the  input  signal  after  being  split  in  I  and  Q  channels 
is  weighted  with  the  optimum  set  of  weights  A  and  B  for  each 
channel.   The  results  are  then  combined  in  order  to  get  an  out- 
put of  positive  doppler  as  well  as  an  output  of  negative  doppler, 
Fig.  37  represents  the  implementation  of  the  filter  that  gener- 
ates the  sets  of  A  and  B  weights  for  each  of  the  M  frequency 
sub- intervals  in  a  channel.   The  T.'s  are  the  interpulse  delays 
and  the 

,        0  <  i  <  N  -  1 
a  .  .  b  .  . 

1J  '  1J      0  <  j  <  M  -  1 

are  the  coefficients  determined  by  the  optimization  process. 
It  should  be  noted  that  in  this  optimization  process  the 
interpulse  periods  are  not  optimized:  only  the  weights  are 
optimized.   Fig.  38  relates  the  signal-to-clutter  gain  aver- 
aged over  the  entire  velocity  region  with  the  signal-to-gain 
maximized  in  each  of  eight  frequency  sub- intervals .   Both 
N  =  5  and  N  =  10  are  presented.   A  substantial  improvement 
due  to  optimization  is  evident.   Another  optimization  tech- 
nique is  developed  in  [18]. 
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The  optimized  processors  are  called  Constrained  Improvement 
MTI  Processors  (CIP) .   Given  a  specified  improvement  factor, 
the  number  of  pulses  to  be  processed  and  the  PRF  stagger 
sequence,  the  mean  square  deviation  from  a  constant  response 
in  the  velocity  region  is  minimized.   Two  main  approaches  can 
be  given  to  the  problem.   In  the  first,  the  PRF  stagger 
sequence  is  fixed  and  the  weighting  coefficients  are  chosen 
in  order  to  optimize  the  processor;  in  the  second,  the  weight- 
ing coefficients  are  fixed  while  the  stagger  sequence  is  chosen 
in  order  to  optimize  the  processor.   Defining 

F  * 
f  =  tt-  ,    F'  =  maximum  PRF 

F   =  unstaggered  PRF 

figures  59,  40  and  41  represent  the  optimum  filter  response 
with  a  four  pulse  return,  f  =  8,  o  -  0.01  and  I  =  30  db , 
using  respectively  a  linear  PRF  stagger  with  ±201  interpulse 
variation,  and  a  sinosoidal  PRF  stagger  with  ±10%  (Fig.  40) 
and  ±90%  (Fig.  41)  interpulse  variation.   As  it  can  be  seen, 
as  the  variation  increases  the  response  becomes  more  uniform. 
Comparing  figures  42,  39  and  43,  in  which  the  number  of  pulses 
are  respectively  3,  4  and  6  for  a  common  interpulse  variation 
of  20%,  it  can  also  be  seen  that  the  response  improves  with 
the  number  of  pulses  processed. 

A  comparison  between  optimization  procedures,  different 
from  a  direct  comparison  between  responses,  can  be  devised 
using  the  fraction  of  frequencies  for  which  a  response  is 
less  than  some  specified  value  as  a  comparison  parameter. 
If  the  0  db  scale  corresponds  to  the  mean  value  of  the 
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Frequency  response  of  four-pulse  CIP  using  sinusoidal  PRF 
Stagger  with  ±10  percent  interpulse  period  variation. 
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frequencies  of  interest,  that  means  that  50%  of  the  frequen- 
cies are  above  as  well  as  below  the  0  db  level.   This  leads 
to  a  cumulative  distribution  function  in  which  instead  of  a 
fraction  of  frequencies,  the  term  probability  is  used.   Fig. 
44  is  a  cumulative  distribution  of  the  responses  of  a  four 
pulse  CIP  and  a  four  pulse  processor  using  the  Ref.  [16] 
optimization  criteria.   This  comparison  criteria  gives  a 
better  performance  for  the  CIP  processor  since  its  curve  is 
closer  to  being  the  ideal  step  response  at  0  db  than  the  other 
one,  but  on  the  average  they  are  similar. 

\   Another  type  of  digital  MTI  processor  is  called  the  matrix 
MTI .   In  this  case  there  is  no  parallel  between  the  analog  and 
the  matrix  MTI  processor.   This  processor  (Fig.  45),  in  addi- 
tion to  the  tremendous  flexibility  in  the  modification  of  the 
shape  of  the  response  curve,  has  associated  with  it  a  very 
simple  clutter  locking  mechanism  that  shifts  the  response 
curve  whenever  there  is  an  average  velocity  associated  with 
clutter.   To  simplify  the  explanation  of  the  processor  in  Fig. 
45,  only  two  levels  of  quantization  are  used.   Basically  the 
processor  is  a  coherent  system  that  compares  the  phase  differ- 
ence between  two  successive  returns  and  weights  that  difference 
in  phase  according  to  the  shape  it  is  intended  to  give  to  the 
response  curve.   Since  d<f>  =  u,  dt,  by  weighting  d4>  in  effect 
the  doppler  is  being  weighted. 

The  phase  angle  of  a  signal  return  is  defined  by 

<J>  =  the"1  9. 
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and  the  amplitude  by  (I2  +  Q2)   .   Since  in  the  case  of 
Fig.  45  there  is  only  one  quantization  level  for  I  and  Q, 
that  means  that  the  only  possible  phases  for  I  are  0°  and 
180°  and  for  Q,  90°  and  270°,  assuming  a  signal  to  noise  ratio 
much  greater  than  1.   When  the  average  phase  difference  is  0°, 
it  implies  a  stationary  target,  but  if  d<$>T   =  180°  that  indicates 
an  optimum  velocity  target.   Since  each  vector  return  can  be 
located  in  all  of  the  four  positions  ±1  ±    jQ,  there  are  4x4 
possible  combinations  for  two  consecutive  returns.   So,  what 
the  compression  matrix  block  in  Fig.  45  does  is  to  generate 
signals  proportional  to  the  phase  difference  output  d(}>T 
between  two  consecutive  vector  returns,  which  in  this  partic- 
ular case  can  take  the  values  0°,  90°,  180°  and  270°.   Each 
of  the  four  outputs  has  a  corresponding  weighting  factor  W 
and  the  result  is  summed  to  generate  the  MTI  output.   Fig.  46 
indicates  the  four  types  with  the  vector  combinations  that 
generated  them,  as  well  as  the  weighting  factors  used  in  order 
to  approximate  an  MTI  response  curve.   The  final  response  curve 
corresponds  to  a  statistical  average  of  d<t>   when  uniformly 
random  phases  are  introduced  at  the  input.   In  order  to  have 
a  signal  proportional  to  the  clutter  velocity,  it  is  only 
necessary  to  subtract  the  outputs  of  the  region  (d)   (Fig.  46) 
from  those  of  region  (c)  and  divide  by  the  number  of  range 
gates,  that  is,  for  two  quantization  levels  and  a  uniformly 

distributed  phase  input 

M  =    [(d)    -     (c)](90°) 

N 

and  VAR   d<{>    =  —    (45°) 

3N 
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So,  in  a  large  clutter  environment  if  N  is  large,  d<fZ   can 

smoothly  adjust  the  local  oscillator  to  place  the  MTI  null 

at  the  mean  value  of  the  clutter  doppler. 

Consider  now  an  n  quantization  bit  circuit.   For  each 

channel  there  are  2   possible  quantization  levels.   If  only 

one  quadrant  of  possible  phase  differences  is  considered,  it 

can  be  seen  that  there  are   (Fig.  47) 

2n    2n    2n 

—  x  —  -  —  +1 

2     2     2 

possible  phases.   The  -2"    +1  parcel  results  from  the  fact 
that  of  all  the  diagonal  combinations  (45°  zone)  only  one 
combination  can  be  counted.   As  an  example,  let's  use  two 
quantization  bits.   That  means  that  there  are 


4[22*-2_  2n-l+  1} 


=  12 
n  =  2 


possible  phases  in  all  four  quadrants,  which  also  implies  12 
possible  phase  differences. 

Figure  48  (a)  represents  the  response  and  figure  48  (b) 
is  the  matrix  output  (phase)  for  each  12  x  12  possible  com- 
bination of  two  consecutive  vector  returns.   Since  there  is 
flexibility  in  modifying  the  weighting  coefficients,  it  is 
possible  to  adapt  the  filter  as  a  function  of  the  type  of 
expected  clutter  spectrum.   Again,  the  mean  clutter  velocity 
can  be  determined  by  subtracting  the  outputs  of  the  (b)  region 
(Fig.  48)  from  those  of  the  (m)  regions  and  applying  the 
result  to  the  VCO  associated  with  the  local  oscillator. 
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B.   DIGITAL  SAR 

The  synthetic  aperture  radar  (SAR)  is  a  typical  example 
of  signal  processing  where  digital  techniques  are  being  sub- 
stituted for  the  former  optical  processors.   With  the  present 
technology  mini- computers  can  compensate  for  the  aircraft 
movement,  and,  as  with  a  real  antenna,  the  synthetic  antenna 
can  be  steered  and  even  scanned.   The  basic  problem  with  SAR 
processors  is  the  need  for  high  storage  capabilities  and  very 
fast  data  rates.   The  film  as  an  optical  storage  device  became 
almost  impractical  since  a  real  time  display  is  impossible. 
The  use  of  storage  tubes  is  also  inadequate  due  to  low  effi- 
ciency and  poor  dynamic  range  and  stability.   With  the  present 
high  speed  A/D  converters  (100  MHZ),  low  cost,  and  high  speed 
compact  digital  storage  devices,  it  becomes  feasible  for  real 
time  SAR  processors. 

Table  IV  shows  a  time  overview  of  the  late  achievements 
in  SAR  techniques. 

Figure  49  is  a  block  diagram  of  a  digital  synthetic 
aperture  radar.   Theoretically,  the  processor  must,  for  each 
range  gate,  perform  the  integral 

/   s(t)  r  (t)  d  t 

AT 

where  s(t)  is  the  signal  return  and  r(t)  the  correlator 
reference  function. 

r(t)  =  AR(t)  exp  -j  (f)R(t) 

AR(t)  =  weighting  function  to  control  the  side  lobes 
of  the  synthetic  antenna  pattern. 

^dC^)  =  phase  reference  that  tracks  the  phase  of  s(t). 
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Date  "*     Development 

1951  Carl-Wiley  postulates  doppler  beam 

sharpening  concept. 

1952  University  of  Illinois  demonstrates 

beam  sharpening  concept. 

1957  First  SAR  imagery  using  optical 

correlator  is  produced. 

Mid  1960's     Analog  electronic  SAR  correlation 

demonstrated  in  non  real  time. 

Late  1960's    Digital  electronic  SAR  correlation 

demonstrated  in  non  real  time. 

Early  1970' s    Realtime  digital  SAR  demonstrated 

with  motion  compensation. 
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Since  <|>R(t)  follows  the  phase  of  s(t)  but  with  an  unknown 
difference,  the  signals  are  processed  in  an  I  (in  phase)  and 
Q  (90°  out  of  phase)  channel  in  order  to  prevent  signal  losses 
The  digital  processor  operates  in  the  following  way:  The  radar 
returns  pass  an  A/D  converter  that  samples  the  signal  at  a 
rate,  at  least  higher  than  the  Nyquist  rate,  and  separates 
the  returns  in  range  bins  which  are  approximately  equal  to 
the  range  resolution.   The  number  of  bits  used  per  range  bin 
are  a  function  of  the  desired  dynamic  range  of  the  processor. 
The  digital  data  passes  through  a  buffer  and  prefilter  which 
translates  the  A/D  rate  to  the  rate  of  the  correlator.   The 
bulk  memory  stores  the  data  corresponding  to  an  integration 
time,  AT.   The  data  is  then  correlated  with  the  reference 
generator  output  in  order  to  form  an  azimuth  line  with  a 
length  corresponding  to  the  number  of  range  gates.   After 
correlation,  new  data  enters  the  memory  and  the  process  re- 
peats itself  to  produce  a  new  azimuth  line. 

From  Fig.  50  it  can  be  seen  that  the  bulk  memory  in  bits 
must  be  equal  to  the  number  of  cells  times  the  average  number 
of  bits  per  cell,  so 

BM  =  bulk  memory  =  2K   N  ND 

a    k 

NR  =  number  of  range  gates  =  — 

r 

A  T 

N   =  number  of  azimuth  data  lines  =  -™-  =   T  f 

The  factor  of  2  is  because  of  the  use  of  an  1  and  Q  channel. 
But  if,  for  example,  the  desired  range  resolution  is  30  m, 
that  corresponds  to  a  bandwidth 
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BW  =  2F"  =  ?  x  ^8  m/gec  =  5.0  MHZ  , 
a    2  x  30  m  ' 

which  implies  an  A/D  sampling  rate  per  channel  of  5.0  MHZ. 
On  the  other  hand,  for  a  1  KHZ  PRF  the  unambiguous  range  is 
80  mi,  which  for  a  typical  AR  =  10  mi  of  mapping,  means  that 
there  is  still  an  equivalent  70  mi  excess  time  to  process  that 
data.   This  results  that  in  fact  there  is  no  need  for  the 
correlator  to  work  at  such  high  data  rates.   Using  a  buffer 
that  accepts  the  data  at  a  very  high  data  rate  during  a  short 
period  of  time  but  delivers  it  at  a  lower  rate  to  the  corre- 
lator, it  is  possible  to  reduce  the  correlator  rates  by  many 
orders  of  magnitude.   In  fact,  analyzing  the  function  of  the 
correlator,  it  can  be  seen  that  it  must  in  AT  seconds  produce 
a  number  of  outputs  equal  to  the  total  number  of  cells,  that 

is,  (L/R  ) ND ,  where  L/r   =  CA  is  the  azimuth  compression  rate, 
a   j\  a 

Since  the  real  antenna  illuminates  (L/r  )Np  cells  and  is  ready 
to  receive  new  data  after  T  seconds,  the  correlator  rate  (CR) 

is 

CR  =  (L/ra^NR   =  CA  ND  f 
j R  r 

defining 

N„  =  K   CA  =  total  number  of  doppler  filters 

Fa  l  l 

K   =  azimuth  over  sample  factor  -  1 


Then 


CR  =  Nc  ND  f 
F   R   r 


if  Nr,  =  500  and  f   =  1  KHZ  for  one  doppler  channel  CR  = 
R  r  ]  F 

.5  MHZ  <  A/D  sampling  rate. 
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On  the  other  hand,  if  the  spacing  corresponding  to  V/T 

is  much  finer  than  the  range  resolution  r  ,  that  means  excessive 

data  is  heing  used.   If  the  PRF  cannot  be  lowered  due  to  power 

or  doppler  considerations,  then  a  prefilter  can  be  used  to 

effectively  reduce  the  PRf  =  f  ,  by  a  factor  f  /f  where  f 
j  r  >      ;  s   r        r 

is  the  prefilter  sampling  rate.   The  prefilter  will  reduce 
the  correlator  rate  by  a  factor  f  /f  .   If  the  required  reduc- 
tion is  in  order  to  have  V/T  of  the  order  of  r  ,  then 

a ' 

c            V  L  1  NF     ......               .  . 

f  ~  —  ~  —  — tf  -  — ?r  >  which  implies  a  new  correlator  rate. 

s    r  r     r  T  '         l 

a  a 

CR  =  K    K   N^  ND  f 
os   s   F   R   r 

K   =  prefilter  over  sample  factor 
os    r  r 

k   =  synthetic  array  weighting 

These  two  factors,  the  bulk  memory  and  the  correlator  rate, 
determine  the  type  of  design  for  a  SAR  processor.   The  main 
objective  is  to  reduce  both  CR  and  BM.   Today's  technology  is 
attacking  this  problem  in  two  ways:  devising  correlator  algor- 
ithms which  will  reduce  the  bulk  memory  and  arithmetic,  at  low 
cost,  with  fewer  power  consuming  elements.   The  first  approach 
is  being  made  by  parallel  and  series  combinations  of  correlator 
channels  and  prefilters,  as  well  as  with  the  use  of  FFT  algor- 
ithms.  Another  method  is  the  use  of  pulse  compression  tech- 
niques but  processed  only  after  the  SAR  processor.   This  will 
reduce  K  which  will  reduce  BM.   The  hardware  improvements 
are  in  CCD  memories  and  LSI  at  much  lower  costs. 

As  it  was  said  before,  with  today's  digital  techniques 
various  mapping  modes  are  possible  as  well  as  motion  compensation 
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Figure  51  shows  the  four  types  of  mapping  used  today.   The 
Squint  mode  (Fig.  51  B)  is  only  a  variation  in  angle  of  the 
side-looking  SAR.   The  doppler  beam  sharpened  mapping  is  a 
PPI  representation  with  an  equivalent  antenna  of  very  high 
azimuth  resolution.   The  spotlight  is  a  mapping  where  a  high 
resolution  snapshot  map  is  generated.   Motion  Compensation 
includes  various  functions:  clutter  tracking,  focusing  accel- 
eration compensation  and  real  antenna  stabilization.   Fig.  52 
is  a  block  diagram  of  a  general  SAR  processor.   The  blocks 
are  basically  the  same  as  those  of  Fig.  49  except  for  the 
motion  compensation  blocks  which  are  now  included. 

To  achieve  lower  data  rates,  various  types  of  algorithms 
are  presently  used.   Lower  rates  will  be  achieved  at  the 
expense  of  more  complex  hardware.   The  general  concept  that 
applies  to  all  of  them  is  the  reduction  of  rates  by  the  in- 
crease in  the  number  of  doppler  filters. 

Figure  53  is  a  block  diagram  of  the  multi-channel  prefilter 

processor  approach  [19]  where  m  is  the  number  of  channels. 

The  prefilter  rate  becomes 

R   -  m  ND  f 
p      R   r 

and  since  the  number  of  filters  Np  is  reduced  byl/m,  the  total 

arithmetic  rate  (TAR)  is 

TAR  =  R  +R  =  a,  m  +  2f 

c    p     1  m 

a  =  Nn  f 
i    R   r 

a  =  K    K   N*  ND/AT 
2    os   s   F   R 

,  d  TAR         a2     n 
and  —j =  a      -    —  =  q 

dm      l  m 
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implies  that  the  value  of  m  that  minimizes  TAR  is 

m     =      rQL2./jV2  =       [K  K         N2/f         Tjl/2 

lcx        os   s   F   r   J 

and  the  minimum  TAR  is 

ot       1/  a       1/  1/ 

TAR  =   a,  (-*-)  /2     +  a    (^)  /2    =    2(a.aJ  '2 


mm 


a  j  2    a  2 


=    2(K        K      f   /AT)1/zNcND 

v  os   s   r    '    F  R 


The  two-stage  correlator  (Fig.  54)  is  another  algorithm  [19]. 
The  idea  is  again  frequency  domain  division  but  only  with  one 
prefilter.  The  equivalent  correlator  rate  will  be  the  sum  of 
the  two  correlator  rates.   So 

TAR  =  (KQs  Ks  NR  NF/AT)  M 

+  K    K   ND  NC/MAT 
os   s   R   F 

where  M  is  the  number  of  first-stage  correlation  channels. 

1/ 
As  in  the  previous  case,  there  is  a  value  of  M  =  (Np)  2  that 

minimizes  TAR.   With  the  FFT  algorithm,  it  is  possible  to 

reduce  even  more  the  total  rates.   Table  V  is  a  summary  of 

the  calculated  values  [19]  of  bulk  memory  as  well  as  total 

arithmetic  rates  for  the  SAR  signal  processing  algorithms. 

Another  approach,  completely  different  from  those  already 

presented  to  solve  the  problem  of  high  data  rates  on  the 

correlator  (Fig.  49)  of  the  SAR  processor,  is  through  the  use 

of  parallel  processing  using  associative  memory  [20].   Fig.  55 

represents  a  block  diagram  of  the  associative  memory  and  the 

related  input/output  registers.   In  the  associative  memory, 
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multi  addressing,  multi  read/write,  multi  logical  and  arith- 
metic operations  are  performed.   The  search/write  and  the 
mask  registers  are  used  to  generate  proper  codes  for  multi 
accessing.   The  operations  are  done  in  a  bit  serial  basis  but 
on  all  words  in  parallel  at  the  same  time.   Ellis  [20]  showed 
that  for  an  X  band  radar  with  range  and  azimuth  resolutions 
of  10  meters,  1000  range  cells,  maximum  range  of  150  km  and 
maximum  aircraft  velocity  of  200  m/sec,  a  50  msec  correlation 
time  is  necessary:  this  implies  a  40  nsec  multiplication  time 
Using  associative  memory,  it  is  possible  [20]  to  perform  the 
operations  of  a  correlation  period  in  31.48  msec.   The  advan- 
tages of  hardware  parallel  operations  are  thus  obvious  over 
the  conventional  process  when  speed  is  an  important  factor. 
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VI.   CONCLUSIONS.  TRENDS 

It  was  shown  that  due  to  improvements  in  integrated 
circuitry,  better  microprocessors,  faster  DFT/FFT  algorithms 
and  lower  cost  of  digital  logic,  most  analog  processors  are 
being  replaced  by  digital  processors.   Very  high  data  rates 
continue  to  be  the  major  problem  of  digital  processors,  but 
with  today's  techniques  of  parallel  processing,  general  pur- 
pose digital  radar  signal  processors  are  already  in  use. 
Despite  the  quantization  noise  inherent  to  digital  processors, 
good  reliability,  extreme  flexibility  and  low  cost  of  digital 
processors  still  give  them  a  tremendous  advantage  over  the 
former  analog  processors. 

The  future  demand  for  more  reliable  and  sophisticated 
radar  systems  will  be  a  function  of  cost  and  military  neces- 
sities.  If  the  cost  of  digital  logic  trends  lower  as  is  its 
present  trend,  and  military  demands  continues  high,  digital 
processors  will  play  an  increasingly  important  role  in  modern 
radar  system's  implementation. 
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